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ABSTRACT 

We present spectroscopic observations and chemical abundances of 16 plane- 
tary nebulae (PNe) in the outer disk of M31. The [O III] A4363 line is detected in 
all objects, allowing a direct measurement of the nebular temperature essential 
for accurate abundance determinations. Our results show that the abundances 



^Partially based on observations obtained with the Apache Point Observatory 3.5-nieter telescope, which 
is owned and operated by the Astrophysical Research Consortium; and on observations obtained at the 
Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science 
Foundation (United States), the Science and Technology Facilities Council (United Kingdom), the National 
Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), Ministerio da 
Ciencia e Tecnologia (Brazil) and Ministerio de Ciencia, Tecnologia e Innovacion Productiva (Argentina). 
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in these M31 PNe display the same correlations and general behaviors as Type II 
PNe in the Milky Way Galaxy. We also calculate photoionization models to 
derive estimates of central star properties. From these we infer that our sample 
PNe, all near the bright-end cutoff of the Planetary Nebula Luminosity Function, 
originated from stars near 2 M©. Finally, under the assumption that these PNe 
are located in M31's disk, we plot the oxygen abundance gradient, which appears 
shallower than the gradient in the Milky Way. 

Subject headings: ISM: abundances, nucleosynthesis, planetary nebulae: general, 
stars: evolution 



1. Introduction 

H II regions and planetary nebulae (PNe) are useful probes of the past and present 
chemical composition of the interstellar medium. The a-element abundances of an ensemble 
of H II regions provide a snapshot of the current element content and radial gradients. 



(Some a-elements are also detectable in B supergiants [e.g. iTrundle et all ( 12002| )]. though 
their galactocentric radius range is more restricted.) On the other hand, leaving carbon 
aside, the a-element abundances found in a planetary nebula reflect those in the ISM at the 
time that the progenitor star formedll] Comparisons of H II and PN abundances at various 
values of galactocentric radius, Rgai, represent the only means of measuring cumulative 
and historical a-element production, providing a glimpse into the chemical histories and 
production processes in the disks of various spiral galaxies. 

Furthermore, plots of a-element abundances against one another constrain the produc- 
tion processes and the relative yields of each element. In this manner the a-element ratios 
can reflect the relative importance of the different a-element enrichment rates by supernovae 
and Wolf-Rayet stars. A comparison of such plots at different values of Rgai may yield indi- 
vidual chemical "fingerprints" for galaxies and the subunits from which they may have been 
assembled. 

In previous papers Henry et al. (2004; HKB04) and Henry et al. (2010; HIO) examined 
these questions for a sample of H II regions and PNe in the solar region of the Milky 
Way (MW). The results showed that halo and disk PNe are distinctive in some of their 
a-element ratios. Also, within the disk, log plots of various a-elements against 0/H show 



^According to stellar evolution theory, the central stars of PNc do not generally enrich the a elements in 
their atmospheres, with minor exceptions to be discussed later. 
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linear trends with slopes near unity, suggesting that although PN abundances may differ 
from one nebula to another, the process(es) that enriched the elements are the same. This 
is profoundly interesting since the progenitor stars of nearby PNe have wide-ranging initial 
masses, compositions, ages and birth sites. 

HIO confirmed that the 0/H ratio shows a measurable slope as a function of Rgai in 
the MW. Importantly, the slopes of log (0/H) vs. Rgai for the PNe and the H II regions 
are consistent if not congruous. HIO also found that the scatter of PNe values from the 
gradient trend is often substantially larger than the measurement uncertainties and the 
scatter in results for H II regions. Estimates of Rgai in the MW are notoriously uncer- 
tain since various distance- measurement methods often disagree, particularly at large Rgai. 
Some of the scatter in the 0/H gradient found by HIO is undoubtedly the result of the 
uncertainties in Rgai. (Interestingly, they also concluded that intrinsic abundance variations 
could contribute significantly to the observed scatter, as we mention in §4.3.) However, the 
distance-determination problem disappears in M31, where all PNe have readily measurable 
values of Rgai, and where the key faint diagnostic emission lines (generally [O HI] A4363) can 
be measure d well in luminous PN e using large telescopes. Propitiously, M31 is rife with PNe; 



a survey by lMerrett et al.l (120061 ) (M06) uncovered 2700 PNe in M31 that are not located in 



streams or other anomalous zones of the galaxy. 



Previous PN abundance studie s in M31 have concentrated on the bulge: e.g.. |Jacoby fc Ciardullo 



(I1999I ) (12 bulge, 3 disk/halo) and [Richer et al.l (Il999[ ) (30 bulge). In contrast, many of the 



PNe identified by M06 lie in the outer regions of the disk where extinction is minimal. Many 
of the PNe seen in projection in the outer disk of M31 are sufficiently luminous that a full 
set of a-element abundances can be measured (with some effort) and compared with those 
of PNe at smaller Rgai and with those in the MW at similar galactocentric distance. These 
PNe in M31 with large and known values of Rgai (and low extinction) provide an opportunity 
to compare the co- dependencies of a-element abundances in this population of PNe to the 
population of PNe near the Sun to see if the processes of a-element abundances are similar; 
this is our primary goal. We note that our abundance results will be directly comparable 
to those for MW PNe since our data calibration and analysis methods are identical to those 
employed in HKB04 and HIO. 

The observations and reductions are summarized in §2. In §3, along with the derived 
nebular abundances, we describe the nebular models and inferred stellar properties. In §4 
we examine abundance correlations, evaluate the population membership of these 16 PNe, 
and discuss the implied oxygen gradient in M31. A summary appears in §5. 
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2. Observations and Reductions 

2.1. Target Selection 

As stated, our primary goal was to measure chemical abundances in M31 PNe and, 
secondarily, to examine the 0/H gradient in M31's disk. Therefore, an initial tactical re- 
quirement of our program was to select PNe that, a) had [OIII] A4363 flux that would be 
detectable using moderately large telescopes; and b) spanned a wide range of galactocentric 
distances. Accordingly we turned to the compilation of PNe in M31 by M06. We selected 
16 bright PNe which are all near the bright-end cutoff of the Planetary Nebula Luminosity 
Function (PNLF). In addition, they are located at sky-projected values of Rgai from 5 to 
about 15 kpc and are not associated with star streams or satellites (Wilkinson, private com- 
munication); they are indicated in Figure 1. Assuming that these PNe lie in the disk (as we 
will discuss in §4.2), this translates to rectified (in-plane) distances between 18 and 43 kpc. 



2.2. Observations 

2.2.1. APO 

Our first step was to observe all targets at Apache Point Observatory (APO). Table 1 
gives details of the observations, along with msoo? and coordinates of the PNe. (All PNe 
lying within its survey footprint are visible in SDSS images in r' and g' filters.) We observed 
these 16 PNe between 2007 and 2010 at APO using the 3.5-m ARC telescope and Dual 
Imaging Spectrograph, which allows simultaneous observation from 3700-9600 A with ~8 
A resolution. All of our spectroscopic data longward of 6400A comes from APO. Since the 
targets are point sources, we observed with slits as well matched as possible to the seeing: 
we used a 6'-long slit that was either 1.5" or 2.0" wide, and we binned on-chip by 2 in the 
spatial dimension for a 1.8 A (blue) or 2.3 A (red) x 0.8" effective pixel size. The observing 
conditions over the APO full- and half-night runs varied from photometric with ~1" seeing 
to variably cloudy with seeing ~ 2.5"; we lost two half-nights to high humidity. APO has no 
atmospheric dispersion corrector, so we attempted to position the spectrograph in real time 
to the parallactic angle in order to minimize such effects. Good fluxes were obtained for all 
of the brighter abundance diagnostic lines for H, He, N, O (except for A4363), Ne, and Ar; 
S was most often detected marginally. 
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2.2.2. Gemini 

The -300 km systemic radial velocity of M31, combined with the bright Hg I A4358 
line from the city of Alamogordo rendered the [0II1]A 4363 line difficult to disentangle at 
APO. The A4363 fluxes for all but two of the target PNe were obtained using Gemini North 
where A4358 is invisible (see below). 

The spectral range from 3600 to 6400 A was reobserved at Gemini North in October 2009 
where we observed 14 of the 16 APO PNe using GMOS. PN6 and PNIO were not included 
in these observations, and all results for them come from the APO observations. GMOS has 
~7A resolution. We used a 5.5'-long slit that was 1.5" wide, binning on-chip by 4 pixels 
in the wavelength dimension and by 2 pixels in the spatial dimension, yielding an effective 
pixel size of 1.8 A X 0.15". We were able to measure A4363, unaffected by the Hg line, in all 
14 PNe; the analogous [N II] line, at A5755 was detected in 13 objects. Interestingly, [O III] 
temperatures derived from these higher signal-to-noise observations were almost uniformly 
lower (weaker A4363) than from APO observations, confirming a well known tendency for 
weakly-detected lines to be systematically over-measured when accurate sky subtraction is 
difficult. Observing conditions at Gemini were excellent: conditions were photometric on all 
three nights, and seeing varied between 0.5" and 1.3." GMOS has no atmospheric dispersion 
corrector, so, as at APO, we positioned the spectrograph to the parallactic angle during 
the night. Compared with M06 our observed values of m (5007)1 are systematically slightly 
fainter and agree to within ±0.07 magnitudes (excluding PN6 and PNll, which were poorly 
centered in the spectrograph slit at APO and Gemini, respectively.) 



2.3. Data Reduction 



Data reduction was carried out in with tasks in the noao and the gemini pack- 

ages for DIS and GMOS data, respectively. The two-dimensional PN and standard star 
frames were bias subtracted, flat-field corrected, and wav elength calibrated. Cosmic rays 
were removed using the task lacosmic (Ivan Dokkumll200ll ). One-dimensional spectra were 
extracted and corrected for atmospheric extinction; flnally, a standard star flux calibration 
was applied to the PN spectra, which were then combined into a flnal spectrum for each 
object. Figure 2 shows an expanded Gemini spectrum of PNl with important lines labelled; 



^m(5007)=-2.51ogF(5007)-13.74 (jjacobv 



198' 



i) 



■^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associ- 
ation of Universities for Research in Astronomy under cooperative agreement with the NSF 
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it is generally representative of all the target spectra. 



Results 



3.1. Line Intensities 



We measured emission-line fluxes from the flnal combined, calibrated blue and red spec- 
tra of each PN using splot. These flu xes formed the input for our abundance determinations 
using ELSA, our 5- level atom code fj Johnson et. alll2006l ). The flrst step in the analysis is 
to generate a table of line intensities that have been corrected for interstellar reddening and 
for contamination of the hydrogen Balmer lines by coincident recombination lines of He^"'". 
We corrected for the effects of reddening using th e function of Savage & Mathis (1979). 
Details of the analysis using ELSA are described in iMilingo et al.l (120 lOf ). Table 2 contains 
the emission-line measurements. Column entries are as follows: the flrst column lists the ion 
and wavelength designation of each line; f(A) gives the value of the reddening function from 
Savage and Mathis (1979), normalized to H/3=0; F(A) is the measured flux and estimated 
error, relative to H/3=100; and 1(A) gives the reddening-corrected intensity and estimated 
error, also relative to H/3=100. At the bottom of each column, for each nebula we list the 
logarithmic reddening parameter, c(H/3), the theoretical Ha/H/? ratio appropriate for the 
nebular temperature and density, and the log of the total observed H/3 flux through the 
spectrograph slit. 



3.2. Plasma Diagnostics and Abundances 

Table 3 contains the plasma diagnostics: electron temperature, T, from [O III] and 
[N II]; and electron density, N, from [S II]. All 16 PNe have T[0 III] calculated from A4363. 
Thirteen have T[N II] calculated from A5755 as well; for the other three PNe, we assume 
T[N II] =10,300 K. In three PNe we were able to measure A6312, the auroral line of [S III], 
so we calculated T[S III] in those cases (though this value is not used in the flnal abundance 
analysis). In only two PNe were both [S II] lines AA6717,6731 detected, enabling a direct 
estimate of the electron density. For most of the remaining objects we detected A6731, 
but not A6717, for which estimated upper limits imply densities in excess of 10,000 cm~^, 
consistent with these being among the brightest PNe in M31. To gauge how the assumed 
density would affect the derived abundances, we carried out calculations for both 10000 cm"-^ 
and 15000 cm~'^. The resulting values differed by less than 10%, except for N/H and N/0, 
which are both ~20% lower at 10000 cm~^. Since the [S II] lines indicate densities above 
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this, we adopt the value of 15000 cm 



Ionic abundances derived using ELSA are given in Table 4. The first column lists the 
ion and wavelength (if more than one line is observed for that ion); the second column 
shows which temperature was used in the calculation; "wm" indicates a weighted mean 
value, (weighted by raw observed flux) for the final ionic abundance to be used in further 
calculation. Ionization correction factors that were derived to compute total abundances 
are lis ted in Table 3; these have been calculated in ELSA as described in iKwitter &: Henry 
( l200ll ). It is important to reiterate that our comparisons of abundances in M31 with those in 
MW PNe (§4.1) are all based on identical methodologies at all stages of the data acquisition, 
calibration and analysis. The to tal elemental abunda nces are sho wn in Table 5; the la st two 
columns give values for the sun (lAsplund et al.ll2009l ) and Orion (lEsteban et al.l 120041 ). 



3.3. Properties of the Central Stars and their Progenitors 

3.3.1. Cloudy Models 



For each PN we calculated a model using version 8.00 of Cloudy (IFerland et al.lll998l ). 
in which we matched a selection of nebular line ratios that are diagnostic of stellar and/or 
nebular properties. Table 6 lists the observed and modeled values of these ratios for each 
PN, along with the resulting parameters, including the luminosity, temperature, and surface 
gravity of the central star (CSPN), and the nebular abundances. We used atmospheres from 



RauchI (120031 ) with log g = 6.5. All models were a spherical shell with an inner radius of 10^'' 
cm; except for two, noted in the table, they were truncated (i.e. matter-bounded) to match 
the observed intensity of the [O II] A3727 line. In all cases, we were able to achieve a good 
match to the observed abundance ratios. Nebular outer radii and implied masses are shown 
in the "Model Parameters" section in Table 6. The radii range between 0.04 and 0.08 pc, 
averaging 0.05 pc, and the implied nebular masses fall between 0.03 and 0.21 Mq, averaging 
~0.07 M0. 

The largest source of uncertainty in the models is the (unknown) abundance of carbon, 
a major nebular coolant. In order to obtain agreement with the observed line intensities 
and nebular temperatures, the required input carbon abundances were many times solar, 
which is probably not realistic. Another ad hoc remedy would have been to assume that 
gas density is a function of radial distance from the central star. But in either case we are 
introducing an unconstrained property into the models. At the same time, we point out 
that the primary goal of the model computations is to determine the central star's effective 
temperature and luminosity. Each of these parameters can be constrained by recombination 
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lines whose strengths are relatively insensitive to electron temperature. For example, the 
ratio of He II A4686 to He I A5876 is positively correlated with stellar effective temperature, 
while the nebular luminosity of either Ha or H/3 tracks the stellar luminosity. To test this 
claim, we chose three PN models in which the C abundance had been raised significantly in 
order to obtain a match with observations (PNe 5, 6, and 11). In each case we reran the 
model after reducing the input C abundance so that the C/0 ratio was equal to the solar 
ratio. The resulting change in He II A4686 to He I A5876 and the luminosity of Ha was 
only a few hundredths of a dex in each case. Therefore, we are confident that in terms of 
stellar parameters, our models remain robust despite the adjustment of an unconstrained 
abundance parameter to help bring the entire set of predicted emission line strengths into 
agreement with those that are observed. 



3.3.2. Analysis of Central Star Properties 



Information regarding the properties of the CSPN and their progenitors can be inferred 
by comparing the positions of CSPNs in a theoretical H-R diagram which also contains 
evolutionary tracks of post AGB stars. From the modeling described above, the bolometric 
luminosity as well as the effective temperature of each CSPN was deduced. This information 
appears in Table 6 and is listed again by object in Table 7. For each object identified in 
the first column of Table 7 we list log Teg and logL/L© in columns 2 and 3. Subsequent 
columns list the fraction of stellar luminosity emitted in A5007, the final (CSPN) mass, Mf, 
the initial mass of the progenitor star. Mi, and the main sequence lifetime, tms- 

Figure 3 is a theoretical H-R diagram, i.e., logL/L p versus logTpff, in which w e show 
five evolutionary tracks from the post- AGB models of IVassiliadis &: Wood I (119941 ). Each 
track is labeled with the mass of the remnant star. In addition, we show the positions of 
our 16 sample objects, based upon the information in the first two columns of Table 7. We 
then estimated the CSPN masses and uncertainties by comparing object positions with the 
model tracks in the figure. The resulting values are those listed in column 4 of Table 7. 

In Figure 4 we plot the extinction constant c(H/3) vs. cor e mass for our objec ts. The 
Milky Way foreground toward M31 imposes E(B-V)=0.062 (jSchlegel et al.l Il998l ). which 
translates into a co ntribution to c(H/3) of 0.088, assuming a ratio of c(H/3)/E(B-V)=1.43 
(IKaler fc Lutzlll985l ). The points in Fig. 4 have had this foreground redd ening removed. The 
figure reveals an evident trend of higher extinction with higher core mass. iJacoby fc CiarduUo 
(119991 ) find similar behavior for their sample of mostly bulge PNe in M31. They note that 
this correlation is expected theoretically, as a larger core mass will evolve from a higher-mass 
progenitor whose evolutionary time scale is more rapid and hence, whose cloud of dust-rich 
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ej ecta still surrounds i t as it becomes a high-excitation PN; this is discussed in more detail 
in lCiarduUo fc Jacobyl (Il999[ ). The linear least-squares slope of the c(H/ j)-core-mass relation 



for ou r data in Fig. 4 is 4.69±0.51, which is similar to values quoted by lCiardullo fc Jacoby 
f ll999h for PNe in the SMC (5.6±0.7) and LMC (6.3±1.3). Their value for the slope in M31 
is 8.5±1.6 (excluding the three most massive cores). We have read off the data carefully from 
the upper left panel in their Fig. 1 (also excluding the three most massive cores), and we 
calculate a slope of 4.75±1.95, about half of their value, and very close to that for our PNe. 
If nothing else, this exercise demonstrates the extreme sensitivity of the calculated slope to 
slight differences in the data values for this small sample. The corrected c(H/3) values for 
our PN sample range fror a 0.01 to 0.36, with a mean of 0.15. This value is lower than the 
0.4 in the LMC quoted by iReid fc PaAe3 J201oh . closer to the 0.28 (i.e., E(B-V)~0.2) found 
bv lHerrmann fc CiarduUol J2009h for M94. and M33 as well as the LMC. 



Figure 5 shows a comparison of the mass distribution of our 16 objects with that of two 
large samples of CSPN and white dwarf masses taken from the literature. The CSPN mass 
distribution for 91 of our MW di s k obj ects,, are shown with thin-line bars, where masses 
were taken from IZhang fc Kwok I (119931 ). Thick-line bars show the spread for a sample of 
247 DA whit e dwar fs with Tg^ > 13, OOOK in the Palomar Green survey as tabulated by 
Liebert et al. ( 2005 1^. Values on the horizontal axis refer to CSPN or white dwarf mass, 
while those on vertical axis indicate the number of objects within each bin, where bin width 
is 0.01 Mq. The CSPN masses of our 16 M31 PNe are indicated with X symbols. (Note 
that the vertical coordinate for these objects is arbitrary and has no meaning). The average 
mass and standard deviation for each sample are indicated. We estimate that the typical 
uncertainty in the M31 CSPN masses is ±O.O4M0 and roughly the same for the other two 
samples. Therefore, values for the average mass for the three groups of objects in Fig. 5 are 
consistent with one another. Interestingly, the MW disk sample comprises about 1/3 nitrogen 
rich Type I PNe, suggesting that the average CSPN mass should be a bit higher than the 
average CSPN in M31. We therefore determined separate averages for the Type I and Type II 
MW disk CSPN and found them to be identical. We conclude that the average values in the 
three samples are indistinguishable. Finally, we wish to point out that although our objects 
were selected based upon bright [O III] A5007 line strengths, this does not necessarily imply 
that the CSPN should have systematically higher masses than those found in the two large 
samples displayed in Fig. 5. As we explain in §4.2 below, L(5007) depends on the location 



iLiebert et all (j2005[ ) claim that DA white dwarfs with Tcfr < 13, OOOK possess convectivc atmospheres 
which bring helium to the surface, resulting in higher inferred masses. Thus, we have followed their lead 
and have ignored all white dwarfs with temperatures below 13, OOOK. None of our own objects falls into this 
category. 
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of the CSPN on its evolutionary track as well as on the age of the nebula. In addition, the 
nebular metallicity and density profile are also likely to play a role in determining [O III] 
luminosity. 

Values for Mi in column 5 of Table 7 were deriv ed using the in i tial-fi nal mass relation: 
Mi = —3.28 + 8.55Mf, where we have taken eq. 1 in [Catalan et al.l ( l2008l ) and solved it for 
Mi. These initial masses range between 1.70 and 2.36 Mq with an average and standard 
deviation of 1.92(±0.19). Column 6 of Table 7 provides an estimate of the progenitors' main 
sequence lifetimes, tms, where these values were determined using model results in Table 45 
of Schaller et al. (1992). The average is 1.34(±.37) Gyr with a range of 0.66-1.8 Gyr. 

Figures 6 and 7 contain plots of seven different abundance ratios, indicated in each panel, 
as functions of stellar mass for our final abundances for the M31 disk objects (Table 5). Error 
bars indicate the abundance uncertainties provided in Table 5 along with uncertainties in 
the progenitor mass given in Table 7. Average values are a lso shown for our samp le (thin 
dot-dashed line) and that of the MW disk (thin dotted line) ( Kwitter &: Henryll2012l ). There 
is close agreement between the abundance averages for M31 and the MW disk in all cases 
except that of N/0. The difference in N/0 averages is no doubt due to the fact that roughly 
1/3 of the PNe in our MW disk sample belong to the nitrogen-rich Type I class, while only 
1 of 16 objects (PN5) in the M31 sample is a Type I. This is consistent with our finding, 
discussed below, that the initial masses implied by the CSPN masses of our objects are below 
the minimum threshold for hot bottom burning and its attendant nitrogen production. 



We also include the model predictions by iKarakas I (120101 . solid bold line) and iMarigo 



(120011 . dashed bold line) in Figs. 6 and 7. (Note that Marigo did not publish predictions of 
neon and sulfur.) The two sets of model predictions for He/H, N/0, and 0/H are in close 
agreement with each other over the stellar mass range relevant to our sample. There is also 
satisfactory agreement between models and observations in the cases of He/H, N/0, and 
Ne/H. 

However, theory and observation are at odds in the cases of 0/H, Ne/H, S/0, Ne/0, 
S/H, and S/ 0. Both Marigo's and Karakas's models were "tuned" to the solar abundance set 
published by lAnders &: Grevesse I (119891 ). In the cases of 0/H, Ne/H, and Ne/0 the theory- 
observation offsets can be accounted for by assuming that the more recent solar abundances 
by Asplund et al. (2009) apply. However, S/H and S/0 offsets cannot be explained in 
this way. Rather, observed sulfur abundances in PNe have been known for several years to 
exhibit levels significantly below expected values, most likely because of a faulty ionization 
correction factor for sulfur (see Henry et al. 2012). Thus, the difference between observation 
and theory for S/H is easily explained by this problem, while that of S/0 is related to 
this same problem along with the updated solar abundances. Therefore, we see no major 
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unexplainable discrepancies between our observed abundances and the predictions of models 
of post AGB stars. 

Finally, our observations indicate that these PNe, all near the bright end of the PNLF, 
share numerous characteristics: in addition to the high [O III] luminosities they were selected 
for, they appear similar in size, density, and some abundance characteristics like N/0 and 
Ne/0. 



4. Discussion 

4.1. Abundance Correlations 

Figures 8-11 display the abundance correlations in our objects along with the corre- 
sponding data for the MW PN samples from HKB04 and HIO. The top plot in Fig. 8 shows 
log(N/0) vs. log(He/H). The M31 PNe are shown as filled circles; MW Type I PNe are 
indicated by open circles, and Type II by open triangles; the sun symbol indicates solar 
values. Note the expected, positive correlation for MW Type I PNe, suggestive of the results 
of hot-bottom burning in more massive progenitors producing nitrogen at the expense of 



oxyge n. This behavior is seen in LMC PNe, but not in SMC PNe (see iKwitter fc Henry 



( 120121 ) for a compilation of abundance correlations in the MW, M31, and the MCs). Nei- 
ther the MW Type II nor the M31 sample exhibits a similar strong trend, arguing for a 
Type II classification for the M31 objects. The only possible exception is PN5, the M31 
PN with the highest N/0, which falls above the majority of the MW Type II objects, and 
into the region occupied by Type I PNe in both this plot and the bottom plot showing log 
N/0 vs. log(0/H). The outlier at low He/H is PNG, (observed only at APO) whose helium 
abundance is low and uncertain. In Figs. 8-11, the correlation between any plotted element 
ratio, X/0, and the measured 0/H value means that if 0/H is underestimated, then X/0 
is overestimated; the opposite is the case if 0/H is overestimated. 

The top plot in Fig. 9 shows log(Ne/H) vs. log(0/H). Here and in the following two 
figures the M31 PNe are shown as solid circles, all MW PNe as smaller open circles, and 
solar values as a sun symbol. The well known close correlation between neon and oxygen 
is evident, supporting the idea that they are produced in the same locations; the slopes 
(~0.97) are identical. In these plots and the following plots for argon and sulfur, the M31 
PNe exhibit a tighter correlation than do the MW PNe; this is likely a reflection of the 
relative homogeneity of the M31 sample compared with MW sample. The lack of an obvi- 
ous correlation between Ne/0 and 0/H, as seen in the bottom fig ure, argues against neon 



production in the progenitors of these PNe (IKwitter fc Henryll2012l ). The behavior of argon 
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shown in the two plots of Figure 10, is similar to that of neon, as expected, albeit with more 
scatter, since the derived argon abundances are less accurate. 

Sulfur is shown in Fig. 11. As evidenced by the MW PNe, the sulfur-oxygen correlation 
is even more diffuse than that of argon-oxygen. In addition, the determination of sulfur 
abundances in the M31 PNe is more challenging than usual: apart from any issues concerning 
the ionization correction factor, the lines from both S"*" and the predominant ionization state, 
S"*"^, are very faint, or absent, in these PNe, and for most of them, the assumed density of 
15,000 cm~^ is a lower limit. Nevertheless, it is clear that the derived sulfur abundances for 
PNe in M31 PNe are below solar, indicating that here, as in the MW and the Magellanic 
Clouds, PNe exhibit the "sulfur anomaly," the systematic underabundance of sulfur relative 
to oxygen (see §3.3.2 and Henry et al. 2012). 



Finall y, we briefly rnention a comparison of some of our results with those of lJacoby fc CiarduUo 



(119991 ) and [Richer et al.l (119991 ). noting again that they con centrated on M31 bulge PNe. The 
N/0 and He/H values found by lJacoby fc CiarduUd (119991 ) extend to higher values than our 
sample, indicating Type I PNe; we found only one possible Type I object. We find that the 
range of oxygen abundances in our sample is very similar to the range found by these au- 
th ors. Also, the relatio i i betw een neon and oxygen for the samples of both, plotted in Fig. 3 
of lJacoby fc CiarduUd (119991 ) overlaps very well with our Fig. 9 (top), again demonstrating 
the near universality of this correlation. 



4.2. Population Membership 



We now explore whether the PNe in the M06 survey are drawn from the outer regions 
of the thin or thick disks of M31. First, the PNe follow the brightness of red starlight along 
the minor axis out to an R-band surface brightness of 24 mag arcsec"^ where most of the 
light arises in the thin disk (M06). Beyond this radius where the thick disk is readily visible 
(Collins et al. 2011, Fig. 9) the statistical scatter in the PN numbers becomes large. The 
dispersions of PN kinematics at large radii, 55 km s~^ (M06), agree at lea st superficially with 
the dispersions of stars in the thick disk, 51 km s~^ (jCoUins et al.l 1201 if ). However, as M06 
point out, a flare in the thin disk like that seen in H I can enhance the velocity dispersion. 

Turning to our target PNe, we noted that t hey are seen at projected R g^i between 5 
and 15 kpc. Using a distance to M31 of 770 kpc (IFreedman fc Madorelll990l ). 1 arcminute 
corresponds to 0.224 kpc. Assuming that the observed PNe are associated with M31's disk, 
we can calculate the rectified Rgai (i.e., distance from the center of M31 in the plane of M31's 
disk) according to: 



- 13 - 



R 



rectified 



(X2 + (Y/cos iffl^, 



where X and Y are the distances along the major and minor axes, respectively (in angu- 
l ar or linear distance ), and % is M31's inclination to the plane of the sky, equal to 77.7° 
(Ide Vaucouleurslll958l ). For each PN we estimated the radial velocity expected for its loca- 
tion if it were following the rotation of the disk (Figure 37 of M06). To do this, we calculated 
6*, the azimuthal angle of the PN's location measured eastward from the major axis in the 
southwest: Q=twL~^ (X/Y). We then calculated 0, the corresponding angle in the plane of 
M31: = 6'/sin %. We next took the difference between the observed and expected radial 
velocities, and divided it by the appropriate value of the velocity dispersion, a (Figure 34 of 
M06). Of these 16 PNe, three have velocities ~3(T from that expected for disk rotation and 
one very near the minor axis differs by 5(j; see Table 8. Given that M06 note substantial 
velocity variations all along the major axis we conclude that these objects are likely members 
of M31's disk population. In that case, their rectified Rgai values range from 18 to 43 kpc. 
For comparison, the window of Rgai for HlO's MW PN sample lies between 5 and 20 kpc. 

Another potential discriminator of PN-stellar-population membership is population age. 
However, the [O III] luminosity criterion that we used to select our targets is likely to select 
only those PNe with the properties that we derived; i.e., final core masses between 0.58 
and 0.7 no matter what the population of stars from which they originate. This is 
shown nicely il l the results of an extensive theoretical study portrayed in Figs. 3 and 4 of 
Mendez et al.l ( l2008l ). Their models show that CSPNe with masses less than ~0.58 Mq 



never form PNe, a result that is empirically confirmed in our Fig. 5 for PNe in the MW and 
M31. Above 0.7 the [O III] peak is so brief (~100y) that the corresponding CSPNe are 
highly unlikely to be represented in any limited sample. Whats more, the fraction of massive 
progenitor stars that form such massive CSPNe is small, and their initial masses so large 
that such stars would dwell in the range of Rgai where H II regions are also plentiful. 

There is additional, albeit circumstantial, evidence that our population of PNe is as- 
sociated with the thin disk. Our PN age estimates can be compared to population ages 
measured from C-M diagrams in the same vicinity. We are unaware of stellar age analyses 
at the outer age of the thin disk where our PNe are sampled. However, the ages of the 
PNe are substantially younger than the ages of the stellar thick disk, 6-8 Gy, or the halo 
( iHammer et al.l (120101 ) and references therein). The mass function of such a population will 
not produce CSPNe above the mass floor of 0.5 8 M^;^. Thus, by inf e rence, our PNe are in the 
outer thin disk where all of the H I is located (IBraun et al.l (120091 ): IChemin et al.l ( 120091 )). 



Their metallicities also suggest an association of our target PNe and the thin disk. 
Fig. 12 indicates that the metallicity of the PNe beyond R25 is only slightly below solar. Ig- 
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noring [0/Fe] gradients, this does not agree particularly well with the metallicity of thin-disk 
stars at about th e same Rrai, [ F e/H] ^ -0.7. The disagreement is even worse for [Fe/H] in the 
thick disk, -1.0 (j Collins et al.l (1201 if ) and references therein). However, t he comparison of 
0/H with Fe/H is complex since [0/Fe] also evolves fIChiappini et al.ll2009l ) at rates that de- 



pend on historical supernova activity, gas infall, and other factors discussed by ICollins et al. 
J201lh . 



Finally, the radial 0/H gradient of PNe matches well with that of H II regions which 
presumably form exclusively in the thin disk (see §4.3 below and Fig. 12). This is probably 
no accident. All considered, it seems much more likely that the PNe formed and evolved in 
the outer parts of the thin disk than in the halo or thick disk. 



4.3. The Oxygen Gradient 



Figure 12 shows the M31 oxygen gradient as exhibited by this sample of PNe; also 
plotted are several M31 population I abundance indicators taken from the literature and 
listed in the legend. For comparison we include PN data for the MW oxygen gradient 
taken from HIO and references therein. HIO derived a slope of -0.058±0.006 dex/kpc, while 
Stanghellini fc HaywoodI f l2010l ) foun d a value of -0.023 ± 0.006. Note that we have plotted 
the abundances as a function of R25 ( iGoodwin et al.lll998l ) to account for differences in disk 
scale length. Using a simple (unweighted) r egression, we find a va lue for the M31 PN gradient 
of -0.011±0.004 dex/kpc. For comparison, iBresohn et al.l ( I2OIOI ) used a weighted regression 
to derive a PN oxygen gradient in M33 of -0.013±0.016 dex/kpc for a sample double the 
size of ours, and IStanghellini et al.l ( I2OIOI ) calculated a linear slope of -0.055±0.02 dex/kpc 
for PNe in M81. 

It is clear from Fig. 12 that the dispersion in 0/H values around the mean M31 PN 
gradient is smaller than the scatter found in local MW PNe. In the MW, distance uncer- 
tainties may account for some of the scatter. We emphasize again that all of the PN data 
shown in Fig. 12 are comparable in quality, the data calibration and reddening corrections 
are all similar, and 0/H abundances were measured from the reddening-corrected line ratios 
in the same way. Among HlO's results regarding the MW oxygen gradient is the conclusion 
that as much as 40% of the scatter about the mean gradient for MW PNe may be intrinsic; 
i.e., due to real differences among PN oxygen ab u ndan ces at the same Rgai, perhaps due 
to hot-bottom burning (HBB); see e.g., iKarakas I ( 120101 ). But this should not be a factor 
for the current sample of M31 PNe, all of whose progenitors have masses below the HBB 
threshold (~ 3 M©). Future observations of PNe spanning an even wider range of Rgai will 
offer a deeper look into the behavior of M3rs oxygen gradient. 
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5. Summary and Conclusions 



We have obtained liigli-quality spectra for 16 PNe in the disk of M31 whose projected 
galactocentric distances range between 5 kpc and 15 kpc, corresponding to rectified distances 
from 18 kpc to 43 kpc. Their spectra contain measurable and clean [O III] A4363 allowing 
direct temperature determinations crucial for accurate abundance determinations. The oxy- 
gen gradient implied by our PN abundances a nd galactocentric distances is similar to that 
found for PNe in M33 bv lBresolin et al.l (120101). but sh allower than the MW PN gradient fit 



of HIO and the M81 result of 



Stanghelhni et all ( 120101 ). 



We have calculated He/H, 0/H, N/0, Ne/0, Ar/0 and S/0 ratios for these PNe and 
compared them with analogous values for MW disk PNe. The M31 PNe are found to display 
the same correlations and general behaviors as Type 11 MW PNe. One exception may be 
PN 5, whose abundance profile and inferred progenitor mass indicate that it is a possible 
Type I. 

Photoionization models using our derived abundances yield central star properties indi- 
cating that these PNe originated from stars with main sequence masses between 1.7 and 2.4 
Mq and ages younger than 2 Gyr. Together with their kinematics, this indicates that they 
belong to M31's extended thin disk. The distribution of core masses in our sa mple appears 



indistinguishable from the distribution of MW PN core masses and from the iLiebert et al. 



(120051 ) white dwarf mass distribution. We examined various abundance ratios as a function 
of progenitor mass and compared them with theoretical predictions. After accounting for the 
effects of updated solar abundances in the models we find reasonable agreement for all cases, 
save S/H and S/0. The PN sulfur deficit has long been intractable, but may be yielding to 
new efforts at correcting it (IHenry et al.ll2012l ). 



The absence of many Type I PNe in our sample is expected since there is no sign of 
the formation of the massive progenitor stars in the ambient stellar population. These PNe 
lie beyond the spiral arms where most massive stars form as evidenced by the paucity of 
H 11 regions with Rgai between 15 and 45 kpc. The computed mass range of the CSPN 
progenitors of 1.70 to 2.36 Mq places them well below the threshold for hot-bottom burning 
and significant nitrogen production during third dredge-up. This is consistent with what we 
see in the observed abundance patt ern, in which 15 of ou r 16 PNe are non-Type I and have 
roughly solar N/0 ratios near 0.17 lAsplund et al.l (|2009[ ). It is also likely that these stars 
will prove to have solar C/0 values, as their mass range is still a bit too low for carbon to 
be significantly enhanced during 3rd dredge-up, according to Karakas's results. 



This work demonstrates the utility of PNe for exploring abundances outside of and 
beyond the spiral arms of galaxy disks, particularly in M31. Bright PNe are not only 
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available in the inner portion of the disk but also out to large galactocentric distances, 
where the density of H II regio ns is low. Though t he identification of H II regions in M31 
has recently been augmented by lAzimlu et al.l (120111 ) whose catalog contains more than 3600 
H II regions extending out to ~24 kpc, these authors confirm that M31 is lacking in the high- 
luminosity H II regions that would be most amenable to abundance studies. Therefore, PNe 
continue to offer the best opportunity for evaluating the chemical evolution of the interstellar 
medium of M31 across its disk. We have explored the spectroscopic properties of a sample 
of PNe near the bright-end cutoff of the PNLF. M06 showed that the PNLF is universal 
throughout M31, so that what we learn from this sample may also apply to similar PNe in 
other parts of M31 and more universally where PNe are forming and evolving. 
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Fig. 1. — M31 PNe observed by Merrett et al. (2006) with our target PNe labeled and 
circled. Dot size indicates relative brightness. Central coordinates of M31 are OOh 42m 
44.3s, +41° 16' 08.5"; horizontal divisions are 5m; vertical divisions are 1°. 
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Fig. 2. — Gemini spectrum (1-hr integration) of PNl, with important hnes labeled. 
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Fig. 3. — Plot of logL/L© versus log Teg, show ing the positions of o u r 16 c entral stars listed 
in Table 8 relative to the model tracks from IVassiliadis &: Wood I ( 1l994l ). where the final 
mass, Mf, corresponding with each track is indicated. 
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Fig. 4. — Extinction constant c vs. core mass for our sample PNe. Foreground Galactic 
extinction has been subtracted off (see text). 
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Fig. 5. — Mass distributio n of 297 whit e dwar f stars in the MW disk from the Palomar 
Green survey presente d by iLiebert et al.l (120051 . thick-hne bars); 91 MW disk CSPN from 



Zhang Sz Kwok I (119931 . thin-hne bars); and CSPN masses of our 16 M31 PNe (each marked 
by an 'X') . Averages and standard deviations are hsted in the legend for each group of 
objects. 
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Fig. 6. — Observed number abundance ratios of He/H, N/0, Ne/0, and S/0 (as indicated 
on the vertical axes) versus progenitor mass, M^, in solar masses for our sample of 16 PNe 
in the disk of M31. Sample averages are given by dot-dashed lines, while the average values 
from our MW disk sample is shown with a dotted l ine. Post-AGB model predictions by 
Marigo ( 2001 . bold dashed lines) and Karakas (j2010 . bold solid lines) are also shown. 
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Fig. 7.— Same as Fig. 6 but for 0/H, Ne/H, and S/H. 
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Fig. 8.— top: log(N/0) vs. log(He/H) for our sample of M31 PNe (filled circles), Milky Way 
disk PNe from HK B04 and HIO (Type I : small open circles; Type II: small open triangles), 
and the Sun from lAsplund et al.l ([20091) (sun symbol); bottom: log(N/0) vs. log(0/H). 
Representative error bars are shown in all plots; in addition, for all log(X/0)-log(0/H) 
plots, the correlation between the two values means that underestimating 0/H increases 
X/0 and overestimating 0/H decreases X/0. 
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Fig. 9.— top: log(Ne/H) vs. Iog(0/H); bottom: log(Ne/0) vs. log(0/H). M31 PNe are 
shown as solid circles, Milky Way disk PNe as small open circles, and the Sun as a solar 
symbol. 
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Fig. 10. — top: log(Ar/H) vs. Iog(0/H); bottom: log(Ar/0) vs. log(0/H). Symbols as in 
Fig. 9. 
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Fig. 11. — top: log(S/H) vs. log(0/H); bottom: log(S/0)vs. log(0/H). Symbols as in Fig. 
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Way (13.4 kpc) were taken from lGoodwin et al.l (119981 ). M31 PNe are shown as filled circles; 
other M31 oxygen probes are shown as well; sources are given in the legend. We include 
Milky Way PNe (small circles) from HKB04, HIO and references therein; the Rgai window for 
these MW PNe lies between 5 and 20 kpc. A representative abundance uncertainty (±0.2 
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- 32 - 



Table 1. Observing Log 



PN 


m(5007)'' 


CMOS Exp (sec) 


DIS Exp (scc)'^ 


RA (2000) 


Dec (2000) 


1-746 


20.45 


3600 


7200 


00:46:18.3 


-1-41:31:09.9 


2-553 


20.54 


4800 


8300 


00:41:54.5 


-1-41:36:16.7 


3-795 


20.56 


3600 


8400 


00:40:26.2 


-1-41:20:27.8 


4-337 


20.58 


2400 


7200 


00:41:41.0 


-1-41:38:14.8 


5-2694 


20.60 


3600 


13200 


00:38:46.2 


-1-41:28:19.6 


6-2471 


20.69 




9955'i''' 


00:43:11.2 


-1-42:20:45.7 


7-2472 


20.69 


3600 


10800 


00:44:02.3 


-1-42:17:16.6 


8-2240 


20.70 


3600 


7200 


00:41:01.8 


-1-40:24:20.4 


9-1596 


20.71 


3600 


7200 


00:40:00.8 


-1-41:04:09.3 


10-1583 


20.71 




15600° 


00:39:22.6 


-1-41:06:57.3 


11-2624 


20.72 


7200 


7800 


00:42:15.1 


-1-41:39:10.1 


12-1985 


20.75 


6600 


3600 


00:37:51.3 


-1-40:45:48.1 


13-2690 


20.77 


7200 


3600 


00:40:58.7 


-1-41:29:32.4 


14-71 


20.83 


7200 


3600 


00:42:59.0 


-1-42:02:14.2 


15-2860 


20.84 


3600 


3600 


00:38:55.0 


-1-41:06:55.3 


16-1074 


20.88 


3600 


3600 


00:40:38.1 


-1-41:16:48.0 



^# preceding the hyphen is an index; # after the hyphen is Merrett et al.'s (2006) ID#; 
''m(5007)= -2.51ogF(5007) - 13.74 

'^DIS exposure in the blue; footnoted if different in the red 
'^DIS exposure in the red: 11155 sec 
"Observed at APO only 
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Table 2a. Fluxes and Intensities 



Line 


f(A) 




PNl 




PN2 




PN3 




PNl 


F(A) 


I{A) 


F{A) 


I{A) 


F(A) 


1(A) 


F(A) 


1(A) 


[O II] A3727 


0.292 


42.3 


45.6±10.15 


35.5 


39.0±8.68 


51.7 


63.6±14.14 


36.1 


44.5±9.88 


He II + H9 A3835 


0.262 










14.9 


18.0±3.82 






[Ne III] A3869 


0.252 


108 


115±24 


100.0 


109±23 


92.8 


111±23 


84.3 


101±21 


He I + H8 A3889 


0.247 


15.4 


16.4±3.41 


13.1 


14.2±2.94 


12.8: 


15.2±5.33: 


10.8 


12.9±2.68 


[Ne III] A3968 


0.225 


34.6" 


36.7±10.40" 


38.3" 


41.2±11.28" 


30.9" 


36.2±10.28" 


33.4" 


39.1±10.85" 


He A3970 


0.224 


14.2" 


15.1" 


13.9" 


15.0" 


12.8" 


15.0" 


12.8" 


15.0" 


He II A4100 


0.188 


0.125" 


0.131" 


3.32(-2)" 


3.53(-2)" 


2.79(-2)" 


3.19(-2)" 






US A4101 


0.188 


23.6" 


24.8±4.72" 


25.1" 


26.7±5.06" 


25.7" 


29.4±5.58" 


20.4 


23.3±4.41 


He II A4339 


0.124 


0.228" 


0.235" 


6.09(-2)" 


6.35(-2)" 


5.24(-2)" 


5.73(-2)" 






H7 A4340 


0.124 


43.2" 


44.6±7.68" 


41.9" 


43.6±7.47" 


42.6" 


46.5±7.97" 


43.3 


47.3±8.09 


[O III] A4363 


0.118 


13.3 


13.7±2.34 


11.8 


12.3±2.09 


6.89 


7.49±1.27 


8.50 


9.25±1.57 


He I A4472 


0.090 


4.26 


4.37±0.71 






4.63:: 


4.93±2.55:: 


4.94:: 


5.26±2.72:: 


He II A4686 


0.036 


12.7 


12.8±1.92 


3.41: 


3.45±1.10: 


3.03 


3.11±0.47 






He I -t- [Ar IV] A4711 


0.030 


3.30 


3.33±0.49 






2.10:: 


2.14±1.10:: 


1.83:: 


1.87±0.96:: 


[Ar IV] A4740 


0.023 


4.27 


4.30±0.63 


2.64:: 


2.66±1.36:: 


2.02:: 


2.06±1.05:: 


1.45::: 


1.48±1.49::: 


He II A4859 


0.000 


0.508" 


0.508" 


0.137" 


0.137" 


0.124" 


0.124" 






up A4861 


0.000 


100" 


100±0" 


100" 


100±0" 


100" 


100±0" 


100 


100±0 


He I A4922 


-0.021 


1.14 


1.14±0.16 










1.44:: 


1.42±0.72:: 


[O III] A4959 


-0.030 


528 


524±71 


480 


476±65 


458 


448±61 


461 


451±61 


[O III] A5007 


-0.042 


1569 


1554±209 


1448 


1429±192 


1351 


1312±176 


1384 


1343±180 


He II A5412 


-0.134 


1.05 


1.02±0.13 














[CI III] A5538 


-0.161 










0.888:: 


0.792±0.400:: 






[N II] A5755 


-0.207 


1.86 


1.77±0.22 


2.14: 


2.00±0.62: 


1.43 


1.23±0.15 


3.06 


2.64±0.33 


He I A5876 


-0.231 


15.3 


14.5±1.80 


18.1 


16.8±2.09 


16.4 


13.9±1.72 


19.2 


16.3±2.02 


[O I] A6300 


-0.313 


4.56 


4.21±0.56 


8.28 


7.48±1.00 


5.93 


4.75±0.63 


11.1 


8.90±1.19 


He II A6311 


-0.315 


4.08(-2)::" 


3.77(-2)::" 






1.15(-2)" 


9.19(-3)" 






[S III] A6312 


-0.315 


2.77::" 


2.56±1.32::" 






0.735" 


0.587±0.080" 


1.59: 


1.27±0.40: 


[O I] A6364 


-0.325 


1.48:: 


1.36±0.69:: 


3.11 


2.80±0.38 


1.92: 


1.53±0.48: 


3.40 


2.70±0.36 


[N II] A6548 


-0.358 


23.1 


21.1±2.97 


23.2 


20.6±2.91 


36.5 


28.3±3.99 


29.4 


22.8±3.22 


He II A6560 


-0.360 


1.76" 


1.61" 


0.488" 


0.434" 


0.506" 


0.392" 






Ha A6563 


-0.360 


313" 


286±0" 


322" 


286±0" 


369" 


286±0" 


369 


286±0 


[N II] A6584 


-0.364 


61.2 


55.9±7.94 


67.1 


59.6±8.47 


83.6 


64.6±9.18 


90.6 


70.0±9.95 


He I A6678 


-0.380 


4.00: 


3.64±1.16: 


3.90 


3.44±0.50 






6.47: 


4.94±1.57: 


[S II] A6731 


-0.389 






5.59: 


4.92±1.57: 






2.36: 


1.79±0.57: 


He I A7065 


-0.443 


8.71 


7.80±1.25 


10.8 


9.36±1.50 


9.50: 


6.94±2.26: 


16.2 


11.8±1.90 


[Ar III] A7136 


-0.453 


10.1 


9.03±1.47 


19.3 


16.6±2.71 


23.3: 


16.9±5.52: 


28.5 


20.7±3.36 


[O II] A7324 


-0.481 






26.3:: 


22.5±11.66:: 






28.4:: 


20.2±10.46:: 


[S III] A9069 


-0.670 






19.6: 


15.7±5.71: 






13.4:: 


8.36±4.52:: 


[S III] A9532 


-0.632 






75.0: 


61.0±21.70: 






59.9 


38.3±8.26 


c 




0.11 




0.14 




0.31 




0.31 




Hq/H/3 




2.86 




2.86 




2.86 




2.86 




log Ttifi^ 




-14.88 




-14.86 




-14.86 




-14.91 





^Deblcndcd. 

'^crgs cm ^ s "'" in our extracted spectra 
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Table 2b. Fluxes and Intensities 



Line 






PN5 




PN6 




PN7 




PN8 


f{A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


[O II) A3727 


0.292 


10.8:: 


13.8±7.44:: 


56.2::: 


60.2±61.44::: 


31.7::: 


35.0±35.74::: 


38.5::: 


47.4±48.28::: 


[Ne III] A3869 


0.252 


82.5 


102±21 


105 


111±23 


95.1 


104±22 


108 


129±27 


He I + H8 A3889 


0.247 


10.3: 


12.7±4.44: 






16.9:: 


18.4±9.80:: 






[Ne III] A3968 


0.225 


26.0° 


31.3±9.29'' 


26.1" 


27.6±8.62" 


25.8:" 


27.9±15.23:" 


42.8" 


50.2±13.32" 


He A3970 


0.224 


12.4" 


15.0" 


14.4" 


15.2" 


14.8: " 


16.0: " 


13.7" 


16.1" 


He II A4100 


0.188 






0.274::" 


0.286::" 






0.132" 


0.151" 


Hc5 A4101 


0.188 


30.0 


35.1±6.66 


17.4::" 


18.2±9.73::" 


28.0 


29.9±5.67 


24.4" 


27.8±5.30" 


He II A4339 


0.124 






0.500" 


0.515" 






0.244" 


0.267" 


H7 A4340 


0.124 


45.0 


49.9±8.54 


46.0" 


47.5±8.22" 


41.4 


43.3±7.42 


45.0" 


49.1±8.45" 


[O III] A4363 


0.118 


10.5 


11.6±1.97 


23.6 


24.3±4.14 


10.8 


11.2±l.gi 


14.4 


15.7±2.67 


He I A4472 


0.090 


10.5: 


11.4±3.70: 


2.13: 


2.18±0.71: 










He II A4686 


0.036 






27.7 


28.0±4.23 






10.9 


11.2±1.68 


He I + [Ar IV] A4711 


0.030 


10.5:: 


10.8±5.53:: 










2.44:: 


2.49±1.28:: 


[Ar IV] A4740 


0.023 


2.47 


2.52±0.37 










3.49:: 


3.55±1.81:: 


He II A4859 


0.000 






1.11" 


1.11" 






0.552" 


0.552" 


H/3 A4861 


0.000 


100 


100 


100" 


100±0" 


100 


100 


100" 


lOO.OiO.OO" 


[O III] A4959 


-0.030 


412 


402±55 


494 


492±67 


432 


427±58 


550 


539±73 


[O III] A5007 


-0.042 


1226 


1185±159 


1523 


1511±204 


1274 


1256±168 


1607 


1560±210 


[N II] A5755 


-0.207 


2.18 


1.83±0.23 






1.83 


1.71±0.21 


2.76 


2.38±0.29 


He I A5876 


-0.231 


18.8 


15.5±1.92 


6.12 


5.81±1.42 


15.4 


14.2±1.76 


17.6 


14.9±1.85 


[O I] A6300 


-0.313 


5.83 


4.49±0.60 






8.35 


7.48±1.00 


9.43 


7.54±1.01 


[O I] A6364 


-0.325 


1.76 


1.34±0.18 






3.35 


2.99±0.41 


3.45: 


2.74±0.86: 


[N II] A6548 


-0.358 


10.5 


7.76±1.09 


17.9 


16.5±2.32 


19.9 


17.5±2.49 


34.9 


27.0±3.80 


He II A6560 


-0.360 






3.82" 


3.52" 






1.96" 


1.52" 


Ha A6563 


-0.360 


386 


286±0 


310" 


286±0" 


321 


283±1 


370" 


286±0" 


[N II] A6584 


-0.364 


52.8 


38.9±5.53 


60.1 


55.4±7.87 


64.7 


57.0±8.15 


104 


80.6±11.45 


He I A6678 


-0.380 


5.77:: 


4.20±2.15:: 






2.66 


2.33±0.34 


4.45:: 


3.39±1.73:: 


[S II] A6716 


-0.387 










2.32 


2.03±0.30 






[S II] A6731 


-0.389 










4.60: 


4.01±1.27: 


4.30:: 


3.26±1.67:: 


He I A7065 


-0.443 


13.4 


9.24±1.48 


5.84:: 


5.29±2.72:: 


9.35 


8.01±1.29 


14.1 


10.3±1.64 


[Ar III] A7136 


-0.453 


16.7 


11.4±1.86 


15.2: 


13.7±4.47: 


15.6 


13.3±2.18 


29.8 


21.6±3.51 


[O II] A7324 


-0.481 










13.9:: 


11.8±6.11:: 






[Ar III] A7751 


-0.539 














14.0:: 


9.52±4.99:: 


[S III] A9069 


-0.670 


30.1:: 


17.2±9.29:: 


25.1: 


21.5±7.81: 


24.5 


19.4±4.47 


49.9: 


30.9±11.23: 


[S III] A9532 


-0.632 


66.1:: 


39.0±20.87:: 






29.0: 


23.2±8.30: 






c 




0.36 




0.10 




0.15 




0.31 




Ha/H^ 




2.86 




2.86 




2.83 




2.86 




log Fh^*" 




-14.87 




-15.19 




-14.94 




-14.98 





^Dcblcndcd. 

'^ergs cm~"^ s""'" in our extracted spectra 
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Table 2c. Fluxes and Intensities 



Line 


f(A) 




PN9 




PNIO 




PNll 




PN12 


F(A) 


I{A) 


F(A) 


1(A) 


F{A) 


1(A) 


F(A) 


1(A) 


[O II] A3727 


0.292 


56.5 


60.7±13.51 


73.3: 


91.2±32.96: 


41.4::: 


47.9±48.81::: 


32.1::: 


42.3±43.14::: 


Ho II + H9 X3S35 


0.262 


18.3 


19.5±4.15 










8.07 


10.3±2.19 


[No III] A3S69 


0.252 


99.5 


106±22 


148 


178±38 


113: 


128±45: 


57.9 


73.5±15.39 


He I + H8 A3889 


0.247 


16.2: 


17.2±6.04: 






29.5 


33.4±7.01 


10.7 


13.6±2.81 


[Ne III] A3968 


0.225 


22. 1'' 


23.4±7.71" 


56.5" 


66.8±16.84" 


27.4" 


30.6±9.29" 


24.5" 


30.3±9.07" 


He A3970 


0.224 


14. 2=^ 


15.0" 


13.8" 


16.4" 


13.7" 


15.3" 


12.1" 


15.0" 


He II A4100 


0.188 


7.37{-2)=' 


7.72(-2)" 


0.389" 


0.448" 


0.472" 


0.519" 






H(5 A4101 


0.188 


20.4" 


21.4±4.07" 


22.1" 


25.4±4.92" 


23.8" 


26.2±5.07" 


24.6 


29.4±5.57 


He II A4339 


0.124 


0.135" 


0.139" 


0.721" 


0.791" 


0.876" 


0.934" 






H7 A4340 


0.124 


41.9" 


43.3±7.43" 


48.6" 


53.4±9.30" 


44.3" 


47.3±8.27" 


43.5 


49.0±8.38 


[O III] A4363 


0.118 


12.5 


12.9±2.19 


24.6:: 


26.8±14.04:: 


14.1 


15.0±2.58 


7.13 


7.97±1.35 


He I A4472 


0.090 














5.17: 


5.63±1.83: 


N III + II A4640 


0.048 










6.62 


6.80±1.04 






Ho II A4686 


0.036 


7.47 


7.54±1.13 


31.8 


32.6±4.95 


49.7 


50.8±7.75 






He I + [Ar IV] A4711 


0.030 










5.24:: 


5.34±2.74:: 


1.15:: 


1.18±0.60:: 


[Ar IV] A4740 


0.023 


3.12:: 


3.14±1.60:: 






9.03:: 


9.17±4.69:: 


0.690:: 


0.705±0.360:: 


He II A4859 


0.000 


0.300" 


0.300" 


1.63" 


1.63" 


2.01" 


2.02" 






up A4861 


0.000 


100" 


100±0" 


100" 


100±0" 


100" 


100±0" 


100 


100 


He I A4922 


-0.021 














1.07 


1.05±0.14 


[O III] A4959 


-0.030 


504 


500±68 


583 


570±78 


554 


548±76 


304 


296±40 


[O III] A5007 


-0.042 


1502 


1488±200 


1908 


1850±250 


1638 


1612±219 


900 


865±116 


He II A5412 


-0.134 










3.97:: 


3.74±1.89:: 






[N II] A5755 


-0.207 


0.890: 


0.847±0.261: 






1.47:: 


1.33±0.68:: 


1.31 


1.07±0.13 


He I A5876 


-0.231 


17.3 


16.3±2.03 


20.0: 


16.8±5.20: 


12.4 


11.2±1.39 


16.9 


13.5±1.68 


[O I] A6300 


-0.313 


4.29 


3.97±0.53 






4.26: 


3.68±1.15: 


4.48 


3.33±0.44 


He II A6311 


-0.315 










0.174:" 


0.150:" 






[S III] A6312 


-0.315 










1.76:" 


1.52±0.52:" 


0.426:: 


0.316±0.160:: 


[O I] A6364 


-0.325 


0.870: 


0.805±0.252: 






1.80: 


1.55±0.48: 


1.76: 


1.29±0.40: 


[N II] A6548 


-0.358 


20.5 


18.8±2.65 


49.9 


38.2±5.43 


19.5 


16.5±2.32 


17.8 


12.7±1.79 


He II A6560 


-0.360 


1.04" 


0.950" 


5.77" 


4.41" 


7.57" 


6.39" 






Ha A6563 


-0.360 


312" 


286±0" 


366" 


280±4" 


339" 


286±0" 


402 


286±0 


[N II] A6584 


-0.364 


64.8 


59.4±8.43 


142 


108±16 


60.8 


51.2±7.27 


49.7 


35.2±5.00 


He I A6678 


-0.380 


4.75:: 


4.33±2.21:: 


2.17: 


1.63±0.52: 






4.65: 


3.25±1.03: 


[S II] A6716 


-0.387 






8.30: 


6.22±1.99: 










[S 11] A6731 


-0.389 


3.70:: 


3.37±1.72:: 


13.8 


10.4±1.54 


6.04 


5.03±0.74 


3.27:: 


2.26±1.16:: 


He I A7065 


-0.443 


5.48: 


4.92±1.60: 










14.2 


9.35±1.50 


[Ar III] A7136 


-0.453 


22.2 


19.9±3.24 


33.8 


24.1±3.96 


22.3: 


18.0±5.88: 


14.2 


9.26±1.51 


[O 11] A7324 


-0.481 






19.5:: 


13.6±7.07:: 










[S III] A9069 


-0.670 


41.0 


34.8±7.93 


86.2: 


52.3±19.08: 


53.8:: 


39.1±21.13:: 






P9 A9228 


-0.610 






12.9: 


8.19±2.89: 










[S III] A9532 


-0.632 






126: 


78.5±28.03: 










c 




0.11 




0.32 




0.21 




0.41 




Ha/HP 




2.86 




2.80 




2.86 




2.86 




log Fh/s'' 




-14.99 




-15.07 




-15.27 




-14.82 





^Dcblcndcd. 

'^crgs cm ^ s "'" in our extracted spectra 
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Table 2d. Fluxes and Intensities 



Line 


f(A) 




PN13 




PN14 




PN15 




PN16 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


F(A) 


1(A) 


[O II] A3727 


0.292 


10.8 


12.0±2.67 


54.2 


64.1±14.27 44.2: 


48.9±17.59: 7.45UL 


8.34±3.00 UL 


He II + H9 A3835 


0.262 














9.18:: 


10.2±5.42:: 


[Ne III] A3869 


0.252 


90.0 


98.6±20.65 


95.3 


110±23 


127 


138±29 


102 


112±23 


He I + H8 A3889 


0.247 






24.6 


28.3±5.88 


15.2:: 


16.6±8.82 


: 18.2: 


20.1±7.04: 


[Ne III] A3968 


0.225 


33. e'' 


36.5±10.34" 


32.3" 


36.7±10.39" 38.2" 


41.3±11.32" 43.8" 


47.7±12.59" 


He A3970 


0.224 


13.8" 


15.0" 


13.2" 


15 0" 


13 9" 


15.0" 


13.8" 


15.0" 


He II A4100 


0.188 






5.51(-2)" 


6.14(-2)" 


8.73(-2)" 


9.32(-2)" 


2.24(-2)" 


2.41(-2)" 


Hi5 A4101 


0.188 






25.4" 


28.2±5.36 


28.8" 


30.7±5.84 


28.0" 


30.1±5.71" 


He II A4339 


0.124 


6.25(-2)'' 


6.54(-2)" 


0.103" 


0.110" 


0.161" 


0.168" 


4.13(-2)" 


4.33(-2)" 


H7 A4340 


0.124 


45.3" 


47.4±8.12" 


45.0" 


48.4±8.30 


" 49.8" 


52.0±8.94 


43.8" 


45.9±7.87" 


[O III] A4363 


0.118 


10.5 


10.9±1.85 


11.5 


12.3±2.09 


12.3 


12.8±2.18 


20.4 


21.4±3.63 


He I A4472 


0.090 


4.49 


4.64±0.75 


6.44 


6.79±1.10 






4.32:: 


4.47±2.31:: 


N III + II A4640 


0.048 










4.21:: 


4.29±2.20 






He II A4686 


0.036 


3.50 


3.55±0.53 


5.88 


6.00±0.90 


8.99: 


9.11±2.92 


2.32:: 


2.35±1.21:: 


He I + [Ar IV] A4711 


0.030 










4.25:: 


4.30±2.20 






[Ar IV] A4740 


0.023 










5.36:: 


5.40±2.76 






He II A4S59 


0.000 


0.141" 


0.141" 


0.239" 


0.239" 


0.362" 


0.362" 


9.37(-2)" 


9.37(-2)" 


rip A4oOi 


0.000 


100 


lOOiO 


100 


100±0" 


100" 


100±0" 


100" 


100 itO 


[O III] A4959 


-0.030 


432 


428±58 


442 


434±59 


508 


503±69 


497 


491±67 


[O III] A5007 


-0.042 


1298 


1279±172 


1326 


1295±174 


1526 


1505±202 


1484 


1461 ±196 


[O I] A5577 


-0.170 










4.15:: 


3.91±1.98 






[N II] A5755 


-0.207 






1.88: 


1.67±0.52 


0.776:: 


0.723±0.365:: 0.846:: 


0.782±0.395:: 


He I A5876 


-0.231 


15.8 


14.5±1.80 


17.0 


14.9±1.85 


19.1: 


17.6±5.45 


17.0 


15.6±1.94 


[O I] A6300 


-0.313 


3.84 


3.43±0.46 


8.42 


7.05±0.94 


4.46 


4.01±0.53 


4.69:: 


4.16±2.11:: 


[O I] A6364 


-0.325 






2.26 


1.88±0.25 










[N II] A6548 


-0.358 


1.05: 


0.922±0.291: 


23.4 


19.1±2.69 


17.2 


15.2±2.15 






He II A6560 


-0.360 


0.509" 


0.447" 


0.928" 


0.755" 


1.30" 


1.15" 


0.340" 


0.296" 


Ha A6563 


-0.360 


326" 


286±0" 


351" 


286±0" 


324" 


286±0" 


328" 


286±0" 


[N II] A6584 


-0.364 


26.0 


22.8±3.24 


64.6 


52.5±7.46 


42.2 


37.2±5.29 


30.5 


26.5±3.77 


He I A6678 


-0.380 






6.46: 


5.20±1.65 


4.08:: 


3.58±1.83 






[S II] A6731 


-0.389 


4.35:: 


3.78±1.94:: 


5.36:: 


4.30±2.20 


: 5.64:: 


4.94±2.53 






He I A7065 


-0.443 


7.29:: 


6.21±3.20:: 


6.60: 


5.13±1.67 


8.85: 


7.61±2.47 


7.75:: 


6.54±3.37:: 


[Ar III] A7136 


-0.453 


19.7:: 


16.8±8.65:: 


17.3: 


13.4±4.36 


24.3 


20.8±3.39 


9.05:: 


7.61±3.93:: 


[S III] A9069 


-0.670 






30.3: 


20.7±7.51 


33.2: 


26.4±9.60 






c 




0.16 




0.25 




0.15 




0.17 




Ha/H/3 




2.86 




2.86 




2.86 




2.86 




log Fh^'" 




-14.92 




-14.98 




-15.04 




-15.08 





^Deblended. 

^ergs cm~^ s ^ in our extracted spectra 



Table 3a. Temperatures and Densities 





PNl 


PN2 


PN3 


PN4 




Parameter 


Value Notes 


Value Notes 


Value Notes 


Value 


Notes 



T[oiII] 10900±560 10810±551 9519±427 10060±478 

T[NII] 11910±997 12320±2235 9366±604 13170±1230 

T[oii] ■ ■ • • ■ ■ • ■ ■ 20290±22760 

T[siii] ■ ■ . . ■ ■ • ■ ■ 12000±2162 Used 9532. 

Nojsjj] 15000 Default. 15000 Default. 15000 Default. 15000 Default. 



Table 3b. Temperatures and Densities 







PN5 




PN6 




PN7 




PN8 


Parameter 


Value 


Notes 


Value 


Notes 


Value 


Notes 


Value Notes 


T[OIII] 
T[NII] 
T[oiI] 
Nejsn] 


11290±602 
15220±1660 

15000 


Default. 


13600±873 
10300 

15000 


Default. 
Default. 


11020±648 
12570±4147 
24880±117900 
9193±20440 




11400±613 
11490±926 

15000 Default. 



Table 3c. Temperatures and Densities 





PN9 




PNIO 




PNll 




PN12 




Parameter 


Value 


Notes 


Value 


Notes 


Value 


Notes 


Value 


Notes 


T[oin] 

T[NII] 
T[oiI] 
Ne[s„] 


10830±553 
8356±996 

15000 


Default. 


13270±2768 
10300 

17470±24380 
3854±5037 


Default. 


11090±580 
10790±2702 

15000 


Default. 


11050±576 
11680±958 

15000 


Default. 



Table 3d. Temperatures and Densities 





PN13 




PN14 




PN15 




PN16 




Parameter 


Value 


Notes 


Value 


Notes 


Value 


Notes 


Value 


Notes 


T[oiII] 

T[NII] 

Nc[sil] 


10780±548 

10300 

15000 


Default. 
Default. 


11170±589 

11950±2097 

15000 


Default. 


10760±547 

9438±2045 
15000 


Default. 


13100±810 

11480±3075 

15000 


Default. 
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Table 4a. Ionic Abundances 







PNl 




PN2 




PN3 




PN4 


Ion 




Abundance 




Abundance 




Abundance 


T„sod 


Abundance 


Hc+ 

Hc+2 

icf(Hc) 


[O III] 
fO TTTl 


9.42±1.20(-2) 
1.18±0.18(-2) 
1.00 


[O III] 
[O TTTl 


0.110±0.014 
3. 17itl.02(-3) 
1.00 


[O III] 
fO TTTl 


9.36±1.17(-2) 
2.80±0.43(-3) 
1.00 


[O III] 

[O TTTl 

[w J.J.J.J 


0.109±0.014 
1.00 


O"(6300) 
0^^(6363) 
O" 

0+(3727) 
0+(7325) 
0+ 

0+2(5007) 
0+2(4959) 
0+2(4363) 

icf(O) 


[N II] 
[N III 
wm 
[N II] 
[N II] 
wm 
[O III] 
[O III] 
[O III] 

wm 


*4.92±1.34(-6) 
*4.97±2.79(-6) 
4.93±1.44(-6) 
*2.97±0.79(-5) 

*4.23±1.03(-4) 
*4.12±0.81(-4) 
*4.23±1.03(-4) 
4.20±0.95(-4) 
1.12±0.02 


[N II] 
fN III 
wm 
[N II] 
[N II] 
wm 
[O III] 
[O III] 
[O III] 

wm 


*8.11±3.76(-6) 
* 9.48±4.40(-6) 

8.49±3.89(-6) 
*2.33±1.20(-5) 
*4.43±4.00(-5) 

3.23±2.36(-5) 
*3.99±0.97(-4) 
*3.84±0.75(-4) 
*3.99±0.97(-4) 

3.96±0.S9(-4) 

1.03±0.01 


[N II] 
[N III 
wm 
[N II] 
[N II] 
wm 
[O III] 
[O III] 
[O III] 

wm 


*1.19±0.32(-5) 
* 1.19±0.47(-5) 
1.19±0.32(-5) 
*1.04±0.28(-4) 

*5.61±1.32(-4) 
*5.53±1.04(-4) 
*5.61±1.32(-4) 

5.59±1.21(-4) 

1.03±0.01 


[N II] 
fN III 
wm 
[N II] 
[N II] 
wm 
[O III] 
[O III] 
[O III] 

wm 


*8.81±2.42(-6) 
* 8.36±2.32(-6) 

8.71±2.34(-6) 
*2.36±0.62(-5) 
*3.29±2.14(-5) 

2.77±1.17(-5) 
*4.75±1.13(-4) 
*4.61±0.88(-4) 
*4.75±1.13(-4) 

4.72±1.04(-4) 

1.00 


Ar+2(7135) 
Ar+^(4740) 
icf(Ar) 


[O III] 
[O III] 


*6.99±1.56(-7) 
*5.31±0.94(-7) 
1.19±0.04 


[O III] 
[O III] 


*1.31±0.29(-6) 
*3.37±1.76(-7) 
1.11±0.06 


[O III] 
[O III] 


*1.78±0.63(-6) 
*3.72±1.94(-7) 
1.20±0.06 


[O III] 
[O III] 


*1.91±0.42(-6) 
*2.27±2.30(-7) 
1.06±0.03 


icf(C) 




1.21±0.04 




1.11±0.06 




1.22±0.06 




1.06±0.03 


01+2(5537) 
icf(Ol) 


[O III] 


1.12±0.02 


[O III] 


1.03±0.01 


[O III] 


*1.14±0.60(-7) 
1.03±0.01 


[S III] 


1.00 


N+(6584) 

N+{6548) 

N+(5755) 

N+ 

icf(N) 


[N II] 
[N II] 
[N II] 
wm 


*9.17±2.33(-6) 
*1.02±0.22{-5) 
*9.17±2.33(-6) 
9.44±2.22{-6) 
17.0±4.85 


[N II] 
[N II] 
[N II] 
wm 


*9.38±3.16(-6) 
*9.56±2.98(-6) 
*9.38±3.16(-6) 
9.43±3.06(-6) 
13.7±9.43 


[N II] 
[N II] 
[N II] 
wm 


*1.78±0.45(-5) 
*2.30±0.50(-5) 
*1.78±0.45(-5) 
1.94±0.44(-5) 
6.57±1.73 


[N II] 
[N II] 
[N II] 
wm 


*1.07±0.27(-5) 
*1.03±0.22(-5) 
*1.07±0.27{-5) 
1.06±0.25(-5) 
18.0±7.64 


Ne+2(3869) 
Ne+2(3967) 
icf(Ne) 


[O III] 
[O III] 


*8.45±1.99(-5) 
8.93±2.57(-5) 
1.20±0.04 


[O III] 
[O III] 


*8.23±1.93(-5) 
1.04±0.29(-4) 
1.09±0.04 


[O III] 
[O III] 


*1.39±0.32(-4) 
1.51±0.43(-4) 
1.22±0.06 


[O III] 
[O III] 


*1.01±0.23(-4) 
1.30±0.36(-4) 
1.05±0.02 


S+ 

S+(6731) 
S+2(9532) 
S+2(6312) 
g + 2 

icf(S) 


[N II] 
[N II] 
[O III] 
[O III] 
wm 


*4.34±2.44(-6) 
1.46±0.D9 


[N II] 
[N II] 
[O III] 
[O III] 
wm 


*4.22±1.71(-7) 
4.22±1.71(-7) 
*2.88±l.ll(-6) 

1.50±0.18 


[N II] 
[N II] 
[O III] 
[O III] 
wm 


*1.74±0.44(-6) 
1.24±0.04 


[N II] 
[N II] 
[S III] 
[S III] 
wm 


*1.51±0.53(-7) 
1.51±0.53{-7) 
*1.74±0.57{-6) 
*1.74±0.57{-6) 
1.74±0.57(-6) 
1.54±0.11 
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Table 4b. Ionic Abundances 



PN5 PNG PN7 PN8 



Ion 




Abundance 




Abundance 




Abundance 


Tuacd 


Abundance 


Hc + 


[O III] 


9.92±1.28(-2) 


JO III] 


3.24±0.84(-2) 


[O III] 


9.35±1.40(-2) 


[O III] 


9.50±1.23(-2) 




fO TTTl 




rO TTTl 
[KJ 111] 


9 ^zt+n '^8f 9'i 


fO TTTl 
[W llij 




\Ci TTTl 
[W iiiJ 


1 04 + lf,( 91 


icf{HG) 




1.00 




1.00 




1.00 




1.00 


O"(6300) 


[N II] 


*2.96±0.85{-6) 


[N II] 




[N II] 


*7.59±6.20(-6) 


[N II] 


•9.21zt2.51{-6) 


\^ \_\jO\JO J 


TN" TTl 

[IN iij 


* 9 7fi-l-n RCl( fi\ 
^ . i u HZ u . ou ^ u / 


riV TTl 

[IN lij 




TN" TTl 

[IN llj 


J . zn 1 . 1 rjy ij J 


TN TTl 

[IN iiJ 




o" 




2.91it0.82{-6) 








8.14zb6.62(-6) 




9.54zt2.67(-6) 


0+ (3727) 


fN TTl 

[IN iiJ 


*4.76±2.63(-6) 


fN TTl 


*4.92itl.73(-5) 


fN TTl 

[IN J.J.J 


* 1.46zt3.07(-5) 


fN TTl 

[iM llj 


*3.29±1.28(-5) 


\y y 1 O^O J 


TN" TTl 

[IN iiJ 




TN" TTl 

[IN iiJ 




In TTl 

[IN llj 


* 9 37-1-9 1 8f ^'1 


fN TTl 

[IN iiJ 




0+ 












1.74zt2.66(-5) 






0+^(5007) 


[O III] 


*2.87±0.71(-4) 


[O III] 


*2.06±0.55(-4) 


[O III] 


*3.27ztO.S9(-4) 


[O III] 


*3.68±0.92(-4) 


0+2(4959) 


[O III] 


*2.82±0.56(-4) 


JO III] 


*1.94±0.42(-4) 


[O III] 


*3.22zt0.73(-4) 


[O III] 


*3.67±0.73(-4) 


0+2(4363) 


[O III] 


*2.87±0.71(-4) 


JO III] 


*2.06±0.55(-4) 


[O III] 


*3.27zt0.89(-4) 


[O III] 


*3.68±0.92(-4) 


0+2 


wm 


2.86±0.65(-4) 


wm 


2.03±0.51(-4) 


wm 


3.26zt0.83(-4) 


wm 


3.68±0.84(-4) 


icf(O) 




1.00 




1.79±0.15 




1.00 




1.11±0.02 


Ar+2(7135) 


[O III] 


•8.16±1.84(-7) 


JO III] 


*6.53±2.42(-7) 


[O III] 


*1.01zt0.24(-6) 


[O III] 


*1.51±0.34(-6) 


Ar+^(7751) 


[O III] 




JO III] 




[O III] 




[O III] 


*2.77±1.52(-6) 




wm 




wm 




wm 




wm 


1.91±0.71(-6) 


Ar+^(4740) 


[O III] 


*2.84±0.51(-7) 


JO III] 




[O III] 




[O III] 


*3.90±2.04(-7) 


icf(Ar) 




1.02±0.01 




2.01±0.20 




1.06zt0.08 




1.20±0.05 


icf(C) 




1.02±0.01 




2.22±0.30 




l.OSztO.OS 




1.21±0.05 


icf(Cl) 




1.00 




1.79±0.15 




1.00 




1.11±0.02 


N+{6584) 


[N II] 


*4.55±1.14(-6) 


[N II] 


*9.06±1.58(-6) 


JN II] 


*8.10zt5.72(-6) 


[N II] 


*1.34zt0.34(-5) 


N+{6548) 


[N II] 


*2.66±0.58(-6) 


[N II] 


*7.94±1.38(-6) 


JN II] 


*7.33zt5.10(-6) 


[N II] 


*1.32±0.29(-5) 


N+(5755) 


[N II] 


*4.55±1.14(-6) 


[N II] 




JN II] 


*8.10zt5.72(-6) 


[N II] 


*1.34±0.34(-5) 


N+ 


wm 


4.24±1.05(-6) 


wm 


8.80±1.46(-6) 


wm 


7.93zt5.56(-6) 


wm 


1.33±0.32(-5) 


icf(N) 




61.1±34.61 




9.17±2.66 




19.8zt26.94 




13.5±5.15 


Ne+^(3869) 


[O III] 


*6.54±1.55(-5) 


JO III] 


*3.68±0.92(-5) 


[O III] 


*7.32ztl.86(-5) 


[O III] 


*7.98±1.89(-5) 


Ne+^(3967) 


[O III] 


6.67±2.00(-5) 


JO III] 


3.a3±0.96(-5) 


[O III] 


6.53zt3.61(-5) 


[O III] 


1.05±0.29(-4) 


icf(Ne) 




1.02±0.01 




2.21±0.30 




1.04zt0.09 




1.21±0.05 


S+ 


[N II] 




[N II] 




[N II] 


*2.44±4.73(-7) 


[N II] 


•2.91ztl.56{-7) 


S+(6716) 


[N II] 




[N II] 




JN II] 


2.44zt4.75(-7) 


[N II] 




S+(6731) 


[N II] 




[N II] 




[N II] 


2.44zt4.72(-7) 


[N II] 


2.91±1.56(-7) 


S+^(9532) 


[O III] 


*1.69±0.94{-6) 


[O III] 


*1.63±0.66(-6) 


[O III] 


*2.18zt0.63(-6) 


[O III] 


*3.27ztl.30{-6) 


icf{S) 




2.23±0.58 




1.30±0.06 




1.58zt0.83 




1.39±0.10 
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Table 4c. Ionic Abundances 



PN9 PNIO PNll PN12 



Ion 




Abundance 




Abundance 


Tuacd 


Abundance 


T„,ed 


Abundance 


Hc+ 


rn TTTi 


n 1 n7_i_n m A 


rO TTTl 

[(J iiiJ 


A 1 nt^-l-n O'^ci 
u.xuoznu.uoy 


rn TTTl 
[U iiiJ 


7 91 _(_n Q^f o\ 


m TTTl 


S 7fi-l-1 TOf 91 

0. lOZlZi.iZl^-.ijJ 


He+2 


[O III] 


6.95±1.05(-3) 


[O III] 


2.98±0.46(-2) 


[O III] 


4.68±a.72(-2) 


[O III] 




icf(He) 




1.00 




1.00 




1.00 




1.00 


O*^ (6300) 


[N TTl 

[IN ±±J 


* 1.30±0.58(-5) 


fN TTl 

[iM llj 




fN TTl 

[IN llj 


* 5.42±3.97(-6) 


fN TTl 

[iM iij 


*4.02itl.09(-6) 


UOUO ) 


TN" TTl 

[IN llj 


o . ^ o HZ o ( 


TN" TTl 

[IN iij 




TN TTl 

[IN iij 




fN" TTl 

[IN iij 


* 4 87-1-1 Q9f' fi1 






1.22it0.55{-5) 








5.93it4.16(-6) 




4.26itl.21(-6) 


0''~ (3727) 


fN TTl 


* 1.43it0.76{-4) 


fN TTl 


* 3.88ib2.40(-5) 


fN TTl 


*4.17it3.99{-5) 


fN TTl 

yi\ iij 


*2.88it0.77(-5) 


\J \ 1 OZD ) 


[N II] 




[N II] 




[N II] 




[N II] 




0+ 








4.54ibl.88(-5) 










0+^(5007) 


[O III] 


*4.13±1.01{-4) 


[O III] 


*2.69±1.77(-4) 


[O III] 


*4.15±1.02{-4) 


[O III] 


*2.25±0.55(-4) 


0+^(4959) 


[O III] 


*4.02±0.78{-4) 


[O III] 


*2.39±1.53(-4) 


[O III] 


*4.08±0.80{-4) 


[O III] 


*2.22±0.44(-4) 


0+2{4363) 


[O III] 


*4.13±1.01{-4) 


[O III] 


*2.69±1.77(-4) 


[O III] 


•4.15±1.02{-4) 


[O III] 


*2.25±0.55(-4) 


0+2 


wm 


4.10±0.92{-4) 


wm 


2.62±1.71(-4) 


wm 


4.13±0.94{-4) 


wm 


2.24±0.51(-4) 


icf(O) 




1.07±0.01 




1.28±0.11 




1.65±0.12 




1.00 


Ar+2{7135) 


[O III] 


*1.56±0.35{-6) 


[O III] 


*1.21±0.62(-6) 


[O III] 


•l.34±0.48{-6) 


[O HI] 


*6.95±1.55(-7) 


Ar+^(4740) 


[O III] 


*3.95±2.06{-7) 


[O III] 




[O III] 


*1.08±0.57(-6) 


[O HI] 


*8.40±4.38(-8) 


icf(Ar) 




1.41±0.20 




1.45±0.16 




1.75±0.16 




1.13±0.04 


icf(C) 




1.44±0.21 




1.51±0.20 




1.82±0.22 




1.13±0.04 


icf(Cl) 




1.07±0.01 




1.28±0.11 




1.65±0.12 




1.00 


N+{6584) 


[N 11] 


*1.87±0.65{-5) 


[N II] 


*1.62±0.50(-5) 


[N II] 


*9.64±4.53{-6) 


[N II] 


*5.86±1.50(-6) 


N+{6548) 


[N 11] 


•l.74±0.56{-5) 


[N II] 


*1.68±0.51(-5) 


[N II] 


*9.13±4.12{-6) 


[N II] 


*6.23±1.37(-6) 


N+{5755) 


[N 11] 


*1.87±0.65{-5) 


[N II] 




[N II] 


*9.64±4.53{-6) 


[N II] 


*5.86±1.50(-6) 


N+ 


wm 


1.84±0.62{-5) 


wm 


1.64±0.49(-5) 


wm 


9.52±4.40{-6) 


wm 


5.96±1.42(-6) 


icf(N) 




4.13±1.70 




8.70±3.69 




18.0±15.84 




8.78±2.41 


No+2{3869) 


[O III] 


•7.96±1.87{-5) 


[O III] 


*6.42±4.66(-5) 


[O III] 


*8.80±3.25{-5) 


[O III] 


*5.12±1.20(-5) 


Ne+^{3967) 


[O III] 


5.83±1.89{-5) 


[O III] 


7.99±5.87(-5) 


[O III] 


6.98±2.12{-5) 


[O HI] 


6.99±2.11(-5) 


icf(No) 




1.43±0.21 




1.47±0.22 




1.81±0.22 




1.13±0.04 


S+ 


[N II] 


*5.43±3.17{-7) 


[N II] 


*4.73±2.46(-7) 


[N II] 


*5.04±2.10{-7) 


[N II] 


*2.03±1.09(-7) 


S+(6716) 


[N 11] 




[N II] 


4.77±2.44(-7) 


[N II] 




[N II] 




S+(6731) 


[N 11] 


5.43±3.17{-7) 


[N II] 


4.70±2.48(-7) 


[N II] 


5.04±2.10{-7) 


[N II] 


2.03±1.09(-7) 


S+^(9069) 


[O III] 


•4.05±1.12{-6) 


[O III] 


*4.13±2.38(-6) 


[O III] 


*4.36±2.45{-6) 


[O III] 




S+^(6312) 


[O III] 




[O III] 




[O III] 


*2.40±0.99{-6) 


[O III] 


*5.06±2.81(-7) 


g + 2 


wm 




wm 




wm 


4.30±2.50{-6) 


wm 




icf(S) 




1.17±0.06 




1.32±0.14 




1.48±0.30 




1.29±0.06 
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Table 4d. Ionic Abundances 



Ion 




PN13 




PN14 




PN15 




PNIO 




Abundance 


Tused 


Abundance 


Tuscd 


Abundance 


T„Bcd 


Abundance 


Hc+ 


[O III] 


9.49±1.21(-2) 


[O III] 


9.58±1.23(-2) 


[O III] 


0.115±0.036 


[O III] 


8.99±1.25(-2) 


Hc+^ 


[O III] 


Q T7^n An/ '5^ 
o.z r ±0.4y (^-dj 


[O III] 


0.04ltU.03(^-0j 


[O III] 


43 QC-LO '7A/' 

fS.ooitz. rU(^-o ) 


[O III] 


o T ^ -1-1 1 n/ o\ 
z.lo±i.iU(_-3J 


icf (He) 




1 00 




1 00 




1 00 




1 00 


(6300) 


fN III 


* 5.86±0.86(-6) 


fN III 


*7.97±3.67(-6) 


IN III 


*8.86±6.18(-6) 


fN III 
lit iij 


*4.39±3.72(-6) 


0''(6363) 


[N II] 




[N II] 


*6.64±3.06(-6) 


JN II] 




[N II] 




O" 


wm 




wm 


7.69±3.52(-6) 


wm 




wm 




0+(3727) 


JN II] 


*1.25±0.26(-5) 


[N II] 


*4.05±2.09(-5) 


JN II] 


*7.01±6.18(-5) 


[N II] 


*5.01±4.10(-6) 


0+2(5007) 


[O III] 


*3.61±0.88(-4) 


JO III] 


*3.25±0.80(-4) 


[O III] 


*4.26±1.04(-4) 


[O III] 


*2.23±0.59(-4) 


0+2(4959) 


[O III] 


*3.48±0.68(-4) 


JO III] 


*3.15±0.62(-4) 


[O III] 


*4.12±0.80(-4) 


[O III] 


*2.17±0.46(-4) 


0+2(4363) 


[O III] 


*3.61±0.88(-4) 


JO III] 


*3.25±0.80(-4) 


[O III] 


*4.26±1.04(-4) 


[O III] 


*2.23±0.59(-4) 




wm 


3.58±0.80(-4) 


wm 


3.23±0.74(-4) 


wm 


4.23±0.95(-4) 


wm 


2.22±0.54(-4) 


icf(O) 




1.03±0.01 




1.06±0.01 




1.07±0.03 




1.02±0.01 


Ar+2(7135) 


[O III] 


*1.33±0.72(-6) 


JO III] 


*9.79±3.53(-7) 


[O III] 


*1.66±0.37(-6) 


[O III] 


*3.93±2.14(-7) 


Ar+^(4740) 


[O III] 




JO III] 




[O III] 


*6.91±3.60(-7) 


[O III] 




icf(Ar) 




1.07±0.01 




1.18±0.07 




1.23±0.14 




1.05±0.02 


icf(C) 




1.07±0.01 




1.19±0.07 




1.25±0.17 




1.05±0.02 


icf(Cl) 




1.03±0.01 




1.06±0.01 




1.07±0.03 




1.02±0.01 


N+(6584) 


[N II] 


*4.76±0.75(-6) 


[N II] 


*8.37±2.84(-6) 


JN II] 


*9.15±4.36(-6) 


[N II] 


*3.79±1.79(-6) 


N+(6548) 


[N II] 


*5.67±1.83(-7) 


[N II] 


*8.94±2.81(-6) 


JN II] 


*1.10±0.50(-5) 


[N II] 




N+(5755) 


[N II] 




[N II] 


*8.37±2.84(-6) 


JN II] 


*9.15±4.36(-6) 


[N II] 


*3.79±1.7g(-6) 


N+ 


wm 


4.60±0.74(-6) 


wm 


8.52±2.78(-6) 


wm 


9.68±4.50(-6) 


wm 


3.79±1.79(-6) 


icf(N) 




30.6±7.83 




9.49±4.54 




7.54±5.79 




46.3±37.63 


Ne+2(3869) 


[O III] 


*7.54±1.77(-5) 


[O III] 


*7.35±1.74(-5) 


[O III] 


*1.06±0.25(-4) 


[O III] 


*4.21±1.04(-5) 


Ne+2(3967) 


[O III] 


9.27±2.67(-5) 


JO III] 


8.15±2.35(-5) 


[O III] 


1.06±0.30(-4) 


[O III] 


5.95±1.68(-5) 


icf(Nc) 




1.07±0.01 




1.19±0.07 




1.25±0.16 




1.05±0.02 


S + 


[N II] 


*4.18±2.16(-7) 


[N II] 


*3.72±2.12(-7) 


JN II] 


*6.32±4.13(-7) 


[N II] 




S+(6731) 


[N II] 


4.1S±2.16(-7) 


[N II] 


3.72±2.12(-7) 


JN II] 


6.32±4.13(-7) 


[N II] 




S+2(9069) 


[O III] 




JO III] 


*2.27±0.90(-6) 


[O III] 


*3.12±1.23(-6) 


[O III] 




icf(S) 




1.71±0.14 




1.31±0.10 




1.26±0.14 




1.98±0.64 
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Table 5a. Total Elemental Abundances 



Parameter 


PNl 


PN2 


PN3 


PN4 


Solar Ref 


Orion Ref 


He/H 


0.106±0.012 


0.113±0.014 


9.64±1.17(-2) 


0.109±0.014 


0.085 


0.097 


N/H 


1 61±0 51^41 


1 29±0 65f-4) 


1 27±0 38^4) 


1 91zb0 71(-4) 


6.76(-5) 


5.37(-5) 


N/0 


0.317±0.078 


0.292±0.149 


0.186±0.046 


0.383±0.125 


0.138 


0.100 


0/H 


5.06±1.08(-4) 


4.40±0.97(-4) 


6.83±1.33(-4) 


4.99±1.06(-4) 


4.89(-4) 


5.37(-4) 


Ne/H 


1.02±0.24(-4) 


8.96±2.08(-5) 


1.70±0.38(-4) 


1.06±0.24(-4) 


8.51(-5) 


1.12(-4) 


Ne/0 


0.201±0.032 


0.204±0.034 


0.249±0.039 


0.212±0.035 


0.174 


0.209 


S/H 


6.33±3.67(-6) 


4.95±1.87(-6) 


2.16±0.59(-6) 


2.91±0.98(-6) 


1.32(-5) 


1.66(-5) 


S/0 


1.25±0.67(-2) 


1.13±0.41(-2) 


3.16±0.63(-3) 


5.83±1.96(-3) 


2.70(-2) 


3.09(-2) 


Cl/H 






1.18±0.61(-7) 




3.16(-7) 


2.88(-7) 


Cl/0 






1.73±0.88(-4) 




6.46(-4) 


5.36(-4) 


Ar/H 


1.47±0.23(-6) 


1.83±0.41(-6) 


2.58±0.79(-6) 


2.27±0.52(-6) 


2.51(-6) 


4.17(-6) 


Ar/O 


2.90±0.42(-3) 


4.15±0.84(-3) 


3.78±1.16(-3) 


4.54±0.99(-3) 


5.13(-3) 


7.77(-3) 



Table 5b. Total Elemental Abundances 



Parameter 


PN5 


PN6 


PN7 


PN8 


Solar Ref 


Orion Ref 


He/H 


9.92±1.28(-2) 


5.78±0.63(-2) 


9.35±1.40(-2) 


0.105±0.013 


0.085 


0.097 


N/H 


2.59±1.53(-4) 


8.07±3.07(-5) 


1.57±1.17(-4) 


1.80±0.74(-4) 


6.76(-5) 


5.37(-5) 


N/0 


0.893±0.490 


0.179±0.072 


0.456±0.402 


0.405±0.152 


0.138 


0.100 


0/H 


2.91±0.66(-4) 


4.51±1.03(-4) 


3.44±0.98(-4) 


4.44±0.96(-4) 


4.89(-4) 


5.37(-4) 


Ne/H 


6.65±1.56(-5) 


8.14±2.21(-5) 


7.65±2.24(-5) 


9.64±2.26(-5) 


8.51(-5) 


1.12(-4) 


Ne/0 


0.229±0.036 


0.180±0.029 


0.223±0.036 


0.217±0.034 


0.174 


0.209 


S/H 


3.77±2.41(-6) 


2.12±0.89(-6) 


3.83±1.62(-6) 


4.95±1.94(-6) 


1.32(-5) 


1.66(-5) 


S/0 


1.30±0.80(-2) 


4.70±2.15(-3) 


l.ll±0.60(-2) 


l.ll±0.43(-2) 


2.70(-2) 


3.09(-2) 


Ar/H 


1.12±0.21(-6) 


1.31±0.45(-6) 


1.06±0.27(-6) 


2.76±0.88(-6) 


2.51(-6) 


4.17(-6) 


Ar/O 


3.85±0.65(-3) 


2.91±1.04(-3) 


3.09±0.67(-3) 


6.21±1.96(-3) 


5.13(-3) 


7.77(-3) 
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Table 5c. Total Elemental Abundances 



Parameter 


PN9 


PNIO 


PNll 


PN12 


Solar Ref 


Orion Ref 


He/H 


0.114±0.014 


0.135±0.040 


0.119±0.013 


8.76±1.12(-2) 


0.085 


0.097 


N/H 


7.59±2.28(-5) 


1.42±0.93(-4) 


1.71±1.06(-4) 


5.23±1.66(-5) 


6.76(-5) 


5.37(-5) 


N/0 


0.129±0.043 


0.361±0.125 


0.229±0.153 


0.207±0.051 


0.138 


0.100 


0/H 


5.89±1.34(-4) 


3.94±2.27(-4) 


7.50±1.74(-4) 


2.53±0.53(-4) 


4.89(-4) 


5.37(-4) 


Ne/H 


1.14±0.29(-4) 


9.43±6.28(-5) 


1.60±0.62(-4) 


5.77±1.30(-5) 


8.51(-5) 


1.12(-4) 


Ne/0 


0.194±0.031 


0.239±0.046 


0.213±0.069 


0.228±0.036 


0.174 


0.209 


S/H 


5.37±1.45(-6) 


6.08±3.63(-6) 


7.11±1.72(-6) 


9.15±3.94(-7) 


1.32(-5) 


1.66(-5) 


S/0 


9.12±2.80(-3) 


1.54±0.65(-2) 


9.48±2.35(-3) 


3.62±1.43(-3) 


2.70(-2) 


3.09(-2) 


Ar/H 


2.77±0.63(-6) 


1.76±0.78(-6) 


4.24±1.33(-6) 


8.77±1.77(-7) 


2.51(-6) 


4.17(-6) 


Ar/0 


4.70±0.95(-3) 


4.46±1.13(-3) 


5.65±1.73(-3) 


3.46±0.68(-3) 


5.13(-3) 


7.77(-3) 



Table 5d. Total Elemental Abundances 



Parameter 


PN13 


PN14 


PN15 


PN16 


Solar Ref 


Orion Ref 


He/H 


9.82±1.21(-2) 


0.101±0.012 


0.124±0.036 


9.21±1.26(-2) 


0.085 


0.097 


N/H 


1.41±0.52(-4) 


8.09±2.88(-5) 


7.30±3.48(-5) 


1.76±0.89(-4) 


6.76(-5) 


5.37(-5) 


N/0 


0.367±0.108 


0.210±0.069 


0.138±0.075 


0.757±0.344 


0.138 


0.100 


0/H 


3.83±0.84(-4) 


3.84±0.83(-4) 


5.29±1.26(-4) 


2.32±0.56(-4) 


4.89(-4) 


5.37(-4) 


Ne/H 


8.08±1.88(-5) 


8.75±2.05(-5) 


1.33±0.34(-4) 


4.41±1.08(-5) 


8.51(-5) 


1.12(-4) 


Ne/0 


0.211±0.033 


0.228±0.036 


0.251±0.040 


0.190±0.030 


0.174 


0.209 


S/H 


7.15±2.60(-7) 


3.46±1.28(-6) 


4.73±1.71(-6) 




1.32(-5) 


1.66(-5) 


S/0 


1.87±0.68(-3) 


9.01±3.34(-3) 


8.94±3.53(-3) 




2.70(-2) 


3.09(-2) 


Ar/H 


1.42±0.76(-6) 


1.15±0.41(-6) 


2.88±0.70(-6) 


4.11±2.23(-7) 


2.51(-6) 


4.17(-6) 


Ar/0 


3.72±1.98(-3) 


3.01±1.06(-3) 


5.45±1.18(-3) 


1.77±0.94(-3) 


5.13(-3) 


7.77(-3) 



Table 6a. Observations & Models' 





PNl 


PN2 


PN3 


PN4 


PN5 


PN6 




Obs 


Model 


Obs 


Model 


Obs 


Model 


Obs 


Model 


Obs 


Model 


Obs 


Model 


^OSl-'-A3727+A5007j/-'-H/3 


1 90 

A. .^\J 


1 20 


1 1 7 


1 1 8 


1 14 


1 1 5 


1 14 


1 1 5 


1 08 


1 06 

± .\J\J 


1 20 


1 24 


log Ia3727/Ia5007 


-1.53 


-1.54 


-1.56 


-1.58 


-1.31 


-1.35 


-1.48 


-1.49 


-1.93 


-1.93 


-1.40 


-1.42 


log Ia4363/Ia5007 


-2.05 


-2.05 


-2.07 


-2.04 


-2.24 


-2.24 


-2.16 


-2.11 


-2.01 


-1.98 


-1.79 


-1.81 


log Ia5876/Ih/3 


-0.84 


-0.81 


-0.77 


-0.78 


-0.86 


-0.86 


-0.79 


-0.79 


-0.81 


-0.80 


-1.24 


-1.25 


log Ia4686/Ia5876 


-0.054 


-0.10 


-0.69 


-0.70 


-0.65 


-0.65 




-0.53 




-0.56 


0.68 


0.56 


log Ia6716/Ia6731 




-0.33 




-0.33 




-0.31 




-0.33 




-0.33 




-0.33 


log Ia6584/Ia3727 


0.09 


0.19 


0.18 


0.12 


0.01 


0.06 


0.20 


0.30 


0.45 


0.56 


-0.04 


-0.04 


log Ia3869/Ia5007 


-1.13 


-1.06 


-1.12 


-1.06 


-1.07 


-0.98 


-1.12 


-1.04 


-1.07 


-1.01 


-1.13 


-0.95 


Te (K) 


10900 


10900 


10800 


11100 


9500 


9600 


10000 


10500 


11300 


11600 


13600 


13100 


Model Parameters 


Teff (lO^K) 




120.0 




104.0 




106.5 




109.0 




107.0 




137.0 


log L/Lq 




3.56 




3.60 




3.70 




3.71 




3.85 




3.35 


Ne (cm-3) 




15000 




15000 




15000 




15000 




15000 




15000 


Mneb (Mq) 




0.042 




0.047 




0.071 




0.064 




0.084 




0.026 


R(pc) 




0.039 




0.040 




0.043 




0.042 




0.044 




0.037 


He/H 




0.101 




0.100 




0.085 




0.100 




0.095 




0.056 


0/H (xlO^) 




5.06 




4.40 




6.82 




4.99 




2.91 




3.48 


N/0 




0.32 




0.29 




0.19 




0.38 




0.89 




0.23 


Ne/0 




0.20 




0.20 




0.25 




0.21 




0.23 




0.23 


S/0 (xl02) 




1.25 




2.03 




0.31 




0.57 




1.29 




0.61 


Ar/0 (xlO^) 




2.90 




4.16 




3.78 




4.55 




3.85 




3.76 



"All central star models have log g=6.5; all models except for PN3-795 were truncated to match the observed 13727- 



Table 6b. Observations & Models' 





PN7 


PN8 


PN9 


PNIO 


PNll 


PN12 




Obs 


Model 


Obs 


Model 


Obs 


Model 


Obs 


Model 


Obs 


Model 


Obs 


Model 


log(lA3727+A5007)/lH/3 


1 1 1 


1 1 1 




1 21 


1 19 


1 1 9 


1 29 


1 27 


1 22 


1 22 


96 


97 


T /T 

iOg iA3727/ ■'■A5007 


-1.56 


-1.55 


-1.52 


-1.53 


-1.39 


-1.40 


-1.31 


-1.28 


-1.53 


-1.53 


-1.31 


-1.37 


log Ia4363/Ia5007 


-2.05 


-2.04 


-2.00 


-1.99 


-2.06 


-2.13 


-1.84 


-1.82 


-2.03 


-2.04 


-2.04 


-1.99 


log Ia5876/Ih/3 


-0.85 


-0.82 


-0.83 


-0.81 


-0.79 


-0.78 


-0.78 


-0.77 


-0.95 


-0.90 


-0.87 


-0.87 


log Ia4686/Ia5876 




-0.83 


-0.12 


-0.15 


-0.34 


-0.31 


0.29 


0.28 


0.66 


0.61 




-0.35 


log Ia6716/1a6731 


-0.30 


-0.28 




-0.33 




-0.31 


0.22 


-0.24 




-0.33 




-0.33 


log 1a6584/1a3727 


0.21 


0.10 


0.23 


0.262 


-0.01 


-0.09 


0.17 


-0.11 


0.03 


0.14 


-0.08 


0.01 


log Ia3869/Ia5007 


-1.08 


-1.04 


-1.08 


-1.01 


-1.15 


-1.08 


-1.02 


-0.98 


-1.10 


-1.08 


-1.07 


-0.98 


Te (K) 


11000 


11000 


11400 


11500 


10800 


10400 


13300 


13400 


11100 


11100 


11100 


11500 


Model Parameters 


Teff (lO^K) 




100.0 




119.0 




115.0 




140.0 




182.7 




113.0 


log L/Lq 




3.58 




3.68 




3.46 




3.65 




3.64 




3.84 


Ne (cm-3) 




6300 




15000 




15000 




4300 




15000 




15000 


Mneb (Mq) 




0.116 




0.058 




0.033 




0.211 




0.041 




0.105 


R(pc) 




0.060 




0.041 




0.038 




0.080 




0.039 




0.047 


He/H 




0.093 




0.100 




0.107 




0.129 




0.116 




0.083 


0/H (xlO^) 




3.44 




4.44 




5.90 




3.94 




7.50 




2.53 


N/0 




0.46 




0.41 




0.13 




0.36 




0.23 




0.21 


Ne/0 




0.22 




0.22 




0.19 




0.24 




0.21 




0.23 


S/0 (xl02) 




1.60 




1.333 




0.87 




2.48 




1.19 




0.70 


Ar/0 (xlO^) 




3.08 




6.22 




4.69 




4.48 




5.65 




3.47 



"All central star models have log g=6.5; all models except for PNl 1-2624 were truncated to match the observed 13727- 



Table 6c. Observations & Models^ 





PN13 


PN14 


PN15 


PN16 




Obs 


Model 


Obs 


Model 


Obs 


Model 


Obs 


Model 


iOg(,iA3727+A5007j/iH/3 


1.11 


1.10 


1.13 


1.13 


1.19 


1.21 


1.17 


1.16 


log Ia3727/Ia5007 


-2.02 


-2.02 


-1.31 


-1.31 


-1.49 


-1.49 


-2.24 


-2.24 


log Ia4363/Ia5007 


-2.07 


-2.08 


-2.02 


-2.01 


-2.07 


-2.07 


-2.24 


-2.24 


log Ia5876/Ih/3 


-0.84 


-0.84 


-0.83 


-0.81 


-0.75 


-0.744 


-0.81 


-0.81 


log Ia4686/Ia5876 


-0.61 


-0.59 


-0.40 


-0.39 


-0.29 


-0.34 


-0.82 


-0.78 


log Ia6584/Ia3727 


0.28 


0.05 


-0.09 


-0.02 


-0.12 


-0.18 


0.50 


0.23 


log Ia3869/Ia5007 


-1.11 


-1.07 


-1.07 


-0.97 


-1.04 


-0.96 


-1.12 


-1.07 


Te (K) 


10800 


10700 


11200 


11300 


10800 


10800 


13100 


13400 


Model Parameters 


Teff (lO^K) 




103.5 




112.0 




114.0 




98.0 


log L/Lq 




3.65 




3.58 




3.47 




3.68 


Ne (cm-3) 




15000 




15000 




15000 




15000 


Mneb (Mq) 




0.105 




0.049 




0.032 




0.036 


R(pc) 




0.047 




0.040 




0.038 




0.039 


He/H 




0.091 




0.095 




0.115 




0.089 


0/H (xlO^) 




3.83 




3.84 




5.28 




2.32 


N/0 




0.37 




0.21 




0.14 




0.76 


Ne/0 




0.21 




0.23 




0.25 




0.19 


S/0 (xl02) 




2.61 




0.86 




1.19 




2.72 


Ar/0 (xlO^) 




3.71 




3.00 




5.45 




1.77 



^All central star models have log g=6.5; all models were truncated to match the observed 

13727- 
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Table 7. CSPN and Progenitor Star Properties 



PN # 


logT.ff 


logL/L0 


L5007/L*'' 


Mf/Mo^ 


Mi/MQ<= 


tms(Gyr)d 


1 


5.08 


3.56 


0.13 


0.597 


1.82 


1.53 


2 


5.02 


3.6 


0.12 


0.593 


1.79 


1.62 


3 


5.03 


3.7 


0.13 


0.613 


1.96 


1.17 


4 


5.04 


3.71 


0.11 


0.614 


1.97 


1.12 


5 


5.03 


3.85 


0.09 


0.64 


2.19 


0.85 


6 


5.14 


3.35 


0.09 


0.59 


1.76 


1.67 


7 


5.00 


3.58 


0.09 


0.585 


1.72 


1.76 


8 


5.08 


3.68 


0.12 


0.616 


1.98 


1.12 


9 


5.06 


3.46 


0.12 


0.583 


1.70 


1.8 


10 


5.15 


3.65 


0.12 


0.622 


2.03 


1.07 


11 


5.26 


3.64 


0.10 


0.66 


2.36 


0.66 


12 


5.05 


3.84 


0.05 


0.64 


2.19 


0.89 


13 


5.01 


3.65 


0.08 


0.6 


1.85 


1.48 


14 


5.05 


3.58 


0.11 


0.595 


1.80 


1.57 


15 


5.06 


3.47 


0.12 


0.583 


1.70 


1.8 


16 


4.99 


3.68 


0.06 


0.605 


1.89 


1.39 



^Fraction of modeled central star luminosity observed in A5007 
'^Determined from model tracks of IVassiliadis fc Wood ] {mi). 
Determined from lCatalan et al. I Hooi). 
'^Determined from results of lSchaller et al" ] lll992l '). 
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Table 8. Kinematics of Sample PNe 



PN# 


Rproj (kpc) 


Rrcct (kpc) 


Observed RV (km s"^)"" 


Expected RV (km s"^)"^ 


RV difference in u'^ 


1 


9.6 


23.5 


153 


143 


0.2 


2 


5.0 


19.8 


170 


79 


1.9 


3 


6.0 


23.2 


-33 


-76 


0.9 


4 


5.6 


23.0 


314 


64 


5.2 


5 


10.4 


42.8 


61 


-25^^ 


2.0^^ 


6 


14.8 


37.8 


55 


82*^ 


-0.6^ 


7 


14.1 


29.2 


48 


125^^ 


-1.7^ 


8 


12.4 


20.3 


-86 


-173 


1.8 


9 


7.4 


18.0 


-9 


-162 


3.1 


10 


8.7 


25.1 


-69 


-123 


1.2 


11 


5.3 


18.7 


-33 


113 


-2.9 


12 


14.1 


28.2 


-208 


-132 


-1.6 


13 


5.4 


23.8 


20 


-5 


0.5 


14 


10.3 


27.5 


-20 


119 


-2.9 


15 


9.9 


29.3 


-110 


-104 


-0.1 


16 


5.3 


19.3 


-144 


-105 


-0.8 



'^Value from lMerrett et al.l (|2006l ) corrected for M31's systemic radial velocity of -309 km/sec. 

^Determined from Fig. 37 of iMerrett et al.l ( 20061 ) and corrected for M31's inclination and PN angular 
position. 



^Determined from Fig. 34 of lMerrett et al.l (|2006l ). 
^^Extrapolated from Figs. 34 and 37 of lMerrett et al.l (|2006l ). 



